INGENIEUR-ARCHIV 


UNTER MITWIRKUNG DER 
GESELLSCHAFT FUR ANGEWANDTE MATHEMATIK UND MECHANIK 
ZUSAMMEN MIT 


A. BETZ * K. KLOTTER * K. MAGNUS + E. METTLER 
K. v. SANDEN - E.SCHMIDT - E. SORENSEN 


HERAUSGEGCEBEN VON 


R. GRAMMEL 


XXVI. BAND SECHSTES (SCHLUSS-)HEFT 1958 


SPRINGER-VERLAG - BERLIN / GOTTINGEN / HEIDELBERG 


Abgeschlossen am 15. Dezember 1958 
Postverlagsort Berlin Preis DM 10,60 


INGENIEUR-ARCHIV 


erscheint nach MaBgabe des eingehenden Materials zwanglos in einzeln berechneten Heften, die zu Banden 
vereinigt werden. 
Die fiir das Ingenieur-Archiv bestimmten Manuskripte sind unmittelbar an den Herausgeber 
Herrn Professor Dr.-Ing. Dr. R. Grammel, Stuttgart N, Robert-Bosch-Straf8e 101 
oder an die Herren 
Professor Dr.-Ing. Dr. A. Betz, Gottingen, Herzberger LandstraBe 39A 
Professor Dr.-Ing. K. Klotter, Stanford (Calif.), 621 Alvarado Row. 
Professor Dr. K. Magnus, Stuttgart O, HacklanderstraBe 33 
Professor Dr. E. Mettler, Karlsruhe-Durlach, Geigersberger Str. 12 
Professor Dr.-Ing. K. v. Sanden, Karlsruhe-West, Hertzstr. 16, (T. H. West) 
Professor Dr.-Ing. E. Schmidt, Technische Hochschule, Miinchen, Arcisstr. 21 
Professor Dr.-Ing. E. Sorensen, Augsburg, MAN 
einzusenden. 
Die zum Druck angenommenen Arbeiten werden, soweit dies drucktechnisch méglich ist, nach der 
Reihenfolge ihres Eingangs beim Herausgeber veréffentlicht. 
Die Mitarbeiter erhalten von ihrer Arbeit zusammen 75 Sonderdrucke unentgeltlich. 


Fiir die Abfassung der Arbeiten wird auf das vom Deutschen NormenausschuB herausgegebene Heft 
»Gcestaltung technisch-wissenschaftlicher Veréffentlichungen“ hingewiesen. Die Vorlagen fiir Abbildungen 
sind auf besonderen Blattern erwiinscht und kénnen entweder in Reinzeichnungen oder in klarverstandlichen 
Handskizzen bestehen; die Beschriftung und nétigenfalls die Reinzeichnung nimmt der Verlag vor. 

Mit der Annahme des Manuskriptes erwirbt der Verlag das ausschlieBliche Verlagsrecht fiir alle 
Sprachen und Lander. 


Im ,,Ingenieur-Archiv® erscheinende Arbeiten diirfen vorher an anderer Stelle nicht veréffentlicht 
sein und auch spater nicht anderweitig veréffentlicht werden. 


Photographische Vervielfaltigungen, Mikrofilme, Mikrophote von ganzen Heften, einzelnen Beitragen 
oder Teilen daraus sind ohne ausdriickliche Genehmigung des Verlages nicht gestattet. 


Die Wiedergabe von Gebrauchsnamen, Handelsnamen, Warenbezeichnungen usw. in dieser Zeitschrift 
berechtigt auch ohne besondere Kennzeichnung nicht zu der Annahme, da8 solche Namen im Sinne der Waren- 
zeichen- und Markenschutz-Gesetzgebung als frei zu betrachten waren und daher von jedermann benutzt werden 
diirften, 


SPRINGER-VERLAG 
Heidelberg Berlin- Wilmersdorf 


Neuenheimer Landstrafe 28—30 / Fernsprecher 279 01 Heidelberger Platz 3 / Fernsprecher Sammel-Nr. 83 03 01 


Inhalt: 


Seite 


Paslay, P. R., und A. Slibar, Uber Querschwingungen von gelenkten Anhangern. 
Mit 3 Abbildungen 


Bufler, H., und H.G. Hahn, Kritische Drehzahlen und Biegeschwingungen kon- 
tinuierlich besetzter Wellen bei verschiedenen Lagerungen. Mit 7 Abbildungen . 387 


Zurmiihl, R., Berechnung von Biegeschwingungen abgesetzter Wellen mit Zwischen- 


bedingungen mittels Ubertragungsmatrizen. Mit 4 Abbildungen ............ 398 
Conway, H. D., Some Special Solutions for the Flexural Vibration of Discs of Va- 
rying. Thickness (0653.25.04 sessay (das ee ee eee 408 


Hausenblas, H., Uberschallstrémung am Austritt von Turbinenschaufel-Kranzen. 
Mit 3 Abbildungen 


Trostel, R., Stationdre Warmespannungen mit temperaturabhangigen Stoffwerten. 
Mit 8 Abbildungen |... 0. 7. (.05 5. ses a eetse ise ee 416 


Withum, D., Die Kreiszylinderschale mit kreisférmigem Ausschnitt unter Schub- 
beanspruchung. Mit 8 Abbildungen...............00.0cceccececccceece... 435 


INGENTEUR-ARCHIV 


XXVI.BAND SECHSTES (SCHLUSS-) HEFT 1958 


Uber Querschwingungen von gelenkten Anhangern 
Von P.R. Paslay und A. Slibar 


1. Einfiihrung. Die Querschwingungen von Kraftwagenanhangern wurden in zwei eingehenden 
Arbeiten von H. Ziegler} erstmalig behandelt. Zur Zeit der Durchfithrung der Untersuchungen von 
Ziegler waren keine eindeutigen experimentellen Aussagen iiber die von der StraBenoberflache auf 
einen Luftreifen iibertragenen Krafte bekannt. Von Ziegler wurde daher angenommen, daB die 
auf den Reifen in einer Normalen zur Radebene wirkenden Krafte den Gesetzen der trockenen 
Reibung gehorchen wiirden. H. Ziegler zeigt, da8 — unter Einwirkung der angenommenen 
Reibungskrafte — beim einachsigen Anhanger unter gewissen Verhaltnissen diskrete Werte von 
kritischen Fahrgeschwindigkeiten des Zugfahrzeugs auftreten. Vom schleppenden Fahrzeug wurde 
dabei angenommen, daf es sich langs einer geraden Fahr- 
bahn bewegt, daB aber kleine Querbewegungen des Zug- 
hakens zufolge der Federung des Zugfahrzeugs auftreten 
kénnen. Weiter ergibt die Rechnung, da der zweiachsige 
Anhanger fiir geniigend kleine Fahrgeschwirdigkeiten stabil 
ist und seine Stabilitat gegen Querschwingungen im be- gy 
sonderen durch entsprechende Dampfungsverhaltnisse im 
Zughaken bzw. im Gelenk der Lenkachse verbessert werden — 15 
kann. 

Spater wurde von L. Huber® in grundlegenden Experi- ™ 
menten iiber die von der Strafenoberflache auf die 
Luftreifen itibertragenen Krafte gezeigt, daB die fiir die % 
Fahrtrichtungshaltung und Stabilitat eines Fahrzeugs we- 
sentlichen Seitenfiihrungskrafte nicht nur von der vom 
Reifen aufzunehmenden Vertikallast und einem zugehérigen 
Reibungsbeiwert, sondern wesentlich yom Schraglauf- 
winkel « des luftbereiften Rades abhangig sind. Dieser 
Schraglaufwinkel ist definiert als der vom waagrechten 
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Raddurchmesser und dem Vektor der Relativgeschwindig- Ab. 1. Die Seitenkraft S als Funktion des Schrag- 
- a : laufwinkels « fiir verschiedene Werte der Reifen- 
keit der Radnabe gegen die Strafenoberflache gebildete belastung P. 


Winkel (Abb. 1). M. Olley* hat in einer spiteren Arbeit 
Hubers Ergebnisse bestatigt. Huber und Olley haben experimentell nachgewiesen, daf ein Luftreifen 
nur dann eine Seitenkraft iibertragt, wenn der Vektor der Geschwindigkeit der Radachse nicht in der 
Radebene liegt. Diese von der StraBenoberflache auf den Luftreifen iibertragene Seitenkraft S erwies 
sich als gleichfalls von der aufzunehmenden Vertikallast und der besonderen Ausfiihrungsart eines 
Reifens abhangig. Abb. 1 zeigt charakteristische Werte der Abhangigkeit von S als Funktion des 
Schraglaufwinkels « fiir verschiedene Werte der Normallast P. Mit guter Naherung kann die Un- 
abhangigkeit der Seitenkraft von der Relativgeschwindigkeit des Rades gegeniiber der Strafen- 
oberflache aus den Versuchen gefolgert werden. 
Die experimentellen Ergebnisse von Huber kénnen dargestellt werden durch 
S =~ P.sq(C, 4). (1) 
wenn j einen zwischen Rad und Strafenoberflache auftretenden Reibungsbeiwert, P die vertikale 
Last und C. einen den Reifen kennzeichnenden Druckbeiwert darstellen. Fiir kleine Werte des 
Schraglaufwinkels « kann die Gl. (1) mit Benutzung experimentell bestimmbarer Konstanten ¢ 
und c’ mit guter Naherung auf die Form 
S=(¢—ce P) Pa (2) 
vereinfacht werden. 
1 H. Ziegler, Ing.-Arch. 9 (1938) S. 96 und S. 241. 


2 L. Huber, Dtsche Kraftfahrtforschung 44 (1940) S. 18. 
3 M. Olley, Proc. Inst. Automob. Eng. 41 (1946) S. 147. 


384 P. R. Paslay und A. Slibar, Uber Querschwingungen von gelenkten Anhangern _Ingenieur-Archiv 


Mit Verwendung der Ergebnisse Hubers haben P. Riekert und T. E. Schunck* die rechnerische 
Untersuchung des erforderlichen Lenkwinkels zur vorgeschriebenen Kreisbogenfahrt und der 
resultierenden Bewegung eines Kraftwagens bei vorgegebenem Lenkwinkel untersucht und die 
Stabilitatskriterien angegeben. P. R. Paslay und A. Slibar® haben mit Verwendung der Beziehung 
(1) baw. der Form (2) die rechnerische Bestimmung des Schwimmwinkels, des Lenkwinkels und 
des Leistungsverlustes bei Kurvenfahrt eines Kraftwagens ausgefiihrt und durch Vergleich der 
errechneten Werte mit den von QO. Dietz und R. Harling? durch Versuche gefundenen Werten gute 
Ubereinstimmung von Versuch und Rechnung fiir Lenk- und Schwimmwinkel zeigen kénnen. 

Unter Verwendung der gleichen Gesetze iiber die wirkenden Seitenkrafte an einem Luftreifen 
wurde von den Verfassern? das Verhalten des einachsigen Anhangers untersucht, wenn der Zug- 
haken durch periodisch in das Zugfahrzeug eingefiihrte Steuerbewegungen eine Querschwingung 
ausfiihrt. Es zeigt sich, daB bei Einfiihrung der Seitenkrafte nach den Versuchsergebnissen von 
L. Huber sich keine diskreten Werte fiir bestimmte kritische Fahrgeschwindigkeiten ergeben, 
sondern bei festem Wert von Lenkungsfrequenz des Zugfahrzeugs zu einer charakteristischen 
Frequenz des Anhingers die von den Radseitenkraften gelieferte Dampfung der Vorwartsge- 
schwindigkeit des Anhingers umgekehrt proportional ist. 

2. Die Bewegungsgleichungen. Ein Zugfahrzeug mit Anhanger besitazt viele Freiheitsgrade, 
und eine groBe Anzahl von Parametern wird das Verhalten eines solchen Systems gegen Quer- und 
Vertikalschwingungen beeinflussen. Um die charakteristische Abhangigkeit der Kigenschwingungs- 
zahlen und der Dampfung zu erhalten, wird der Anhanger auf das in Abb. 2 dargestellte System 


Geschwindigkert 
des Zughakens 
\ 


Ersatzrader in Achsmitte 


Abb, 2. Schema des gezogenen Fahrzeugs. 


abgebildet, und es werden zur Aufstellung der Bewegungsgleichungen die Einzelrader durch mittig 
angebrachte Ersatzrader dargestellt. Abb. 3a zeigt die am Modell verwendeten Symbole, Abb. 3b 
das Schema des Lenk- bzw. Lastteiles mit den wirkenden Kraften. Die Behandlung wird fir kleine 
Amplituden der Querbewegung ausgefiihrt und die Zuggeschwindigkeit konstant angesehen.® Da 
die StraBe glatt und eben angenommen wird, kann der Einflu8 der Federung sowohl des Zugfahr- 
zeugs als auch des Anhangers vernachlassigt werden, und der Ersatz der beiden Einzelrader jeder 
Achse durch ein mittig angebrachtes Ersatazrad ist zulassig. Die bei den Querschwingungen auf- 
tretenden Kreiselmomente an den Radern, der Rollwiderstand und das in der Reifenauflage wirkende 

Schraglaufmoment® kénnen bei kleinen Ausschlagen vernachlassigt werden. 

Damit ergeben sich fiir den Lenkteil die Bewegungsgleichungen 
0=X— 5,9 —X, | 
mc,9,= Y—S,—Y, 

ee Loy =—4 Y—4 Xo, —S,(a—e)— Y(h—4)§ —X(L— 4) y,, 

! P. Riekert und T. E. Schunck, Ing.-Arch. 9 (1940) S. 210. 

? P. R. Paslay und A. Slibar, Ing.-Arch. 24 (1956) S. 412. 

° O. Dietz und R. Harling, Dtsche Kraftfahrtforschung 44 (1940). 

' A, Slibar und P. R. Paslay, Journ. Appl. Mech. 24 (1957) S. 515. 

° Mit der Annahme der konstanten Zuggeschwindigkeit wird — um auswertbare Rechenergebnisse zu er- 
halten — bewubt auf die Riickwirkung der Anhangerbewegung auf die Bewegung des Zugfahrzeugs verzichtet. 
Die Beriicksichtigung aller méglichen Freiheitsgrade des Zugfahrzeugs (siehe H. Ziegler, FuBnote 1, S. 383) 
fiihrt bei Benutzung der effektiven Reibungsgesetze des Luftreifens auf ein der Rechnung kaum zugangliches 
Gleichungsystem. Die vorliegende Rechnung wird aber besonders fiir die haufige Kombination eines schweren 
Zugfahrzeugs mit leichtem Anhanger gut brauchbare Ergebnisse liefern. 


® Siehe z. B. E. Fiala, Z. VDI 96 (1954) S. 973. 


(3) 
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wenn X, X, Y, Y die in den Gelenkpunkten auf die Teile des Anhangers iibertragenen Krafte, S, die 

auf das vordere Ersatzrad wirkende Seitenkraft, p, den Querausschlag der Langsachse des Anhanger- 

vorderteils, m, die Masse und I, das Massentragheitsmoment um die Hochachse durch den Schwer- 

punkt des Lenkteils bedeuten. Die Bedeutungen der Langen a, c, und I, gehen aus Abb. 3a hervor. 
Fiir den Lastteil des Anhiangers folgen die Bewegungsgleichungen 


0=X—S. 9, ms (1, Gy + ¢2P2) =Y—S,, I, Pz = — ¢, Y — ¢,X yy — Sy (ly — 9), (4) 
wenn S, die auf das hintere Ersatzrad iibertragene Seitenkraft, ~, den Querausschlag der Langs- 


achse des Anhangerteils, m, und I, dessen Masse bzw. Massentragheitsmoment um die Hochachse 
durch den Schwerpunkt bedeuten. Die Bezeichnungen ¢, und /, sind in Abb. 3a eingezeichnet. 


Abb. 3. a) Geometrie und Geschwindigkeiten am Anhingermodell. 
b) Wirkende Krifte auf Lenk- bzw. Lastteil des Anhingermodells. 


Bezeichnen wir mit v die konstante Fahrgeschwindigkeit des Zugfahrzeugs, so ergeben sich aus der 
Geometrie des Systems nach Abb. 3a fiir kleine Schraglaufwinkel der Ersatzrader die Beziehungen 


age hi tho. +o9 
— mo a ae ia er 2 (5) 


Eliminiert man die im Zughaken T, und im Gelenk T, wirkenden Krafte, so folgt 
(I, + m, cf + mz lt) p, + (m2 ¢g1,)P, +458, +1,8,=0, (6) 
(1, €2 ms) Py + (Lp + mg >) Pp + Sgl, = 0. 
Sind P, und P, die auf den Ersatzradern wirkenden vertikalen Lasten, so liefert die Gl. (2) mit 
Beniitzung von (5) fiir die Seitenkrafte 


S,=(¢—e P,) P,** 


v s v 


nal Fi 7 S. == (c — c 2) Nes I Atl Pe Ar VPs : (7) 


und damit folgen fiir die Beschreibung der resultierenden Querschwingung die beiden simultanen 
Differentialgleichungen 


ee ee es (8) 
ky Py a key Po + kg Qy + hg Pe + hy G2 = 9, 
in welchen die Koeffizienten k,, k,,....k, die Bedeutungen 
bk = Lam, mt, i ted igh Vel ead Aa 
kg => Co Aas ; k, = 1, -+- ms ¢ , | 
pe Pye ee PP, , kg = (c —e' P,) P,, (9) 
Lae ko aC — 6 Pe) Pals 


be (¢-—¢ Pi) P, a, 
haben. 
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Fihrt man in die Bewegungsgleichungen (8) den Lésungsansatz 


eT ey (10 
ein, so folgt 
(ede eas) (aa (11) 
(ka B? + ky B) (hz B? + ke B + hy)/ \B 


und die Sakulargleichung vierter Ordnung zur Bestimmung der Kigenfrequenzen wird 


(ky ky — hi) B* + (hy keg + beg kp — 2 hey hey) B® + (hy hry + Fig keg ar kk, + Keg ke + hi) B? 


+ (hg ky + hts kg — hg hg) B + hg hy = 90. (12) 
Mit den offensichtlichen Bedeutungen von p, qg, r und t kénnen wir (12) in der Form 
Ppp gest rhs oie 2) 


schreiben. 


3. Folgerungen. In (12a) sind q und s die Mafzahlen fiir die Dampfungseigenschaften des 
Anhangers gegen Querschwingungen. Entsprechend dem bekannten Verhalten des einachsigen 
Anhangers! kann erwartet werden, daB auch fiir das System mit zwei Freiheitsgraden der dampfende 
Einflu8 der auf die Luftreifen wirkenden Seitenkrafte mit steigender Fahrgeschwindigkeit des 
Zugfahrzeugs abnehmen wird. Es ist auch aus der Erfahrung bekannt, daB hohe Werte der Fahr- 
geschwindigkeit zu kritischen Fahrzustanden, d. h. zu iibergroBen Queramplituden des Anhangers 
fiihren kénnen?. Fiir die praktische Beurteilung werden daher jene aus (12a) folgenden Eigen- 
kreisfrequenzen von Interesse sein, fiir welche die Dampfung einen kleinen Wert annimmt. Ver- 
nachlassigt man die geringe Verschiebung der Eigenfrequenzen zufolge kleiner Dampfung, so 
miissen beim Auftreten groBer Querschwingungsamplituden die kritischen Lenkungsfrequenzen 
nahe an den ungedampften Eigenschwingungszahlen des Anhangers liegen. Bezeichnet w die 
ungedampften Eigenschwingungszahlen, also die imaginaren Anteile von /, so liefert die Beziehung 


po—reo’®+t=0, (13) 
fiir diese GréBen die Ausdriicke 
if r 1/, 4pt 
O32 = 2p a 2p 1 ae aa a (14) 


wenn die Beiwerte p, r und ¢ in den urspriinglichen Parametern des Systems geschrieben die 
Bedeutungen 


Pp = (A, + m ct + mg lt) (1, + m, ¢3) — (m, c,h)? , 
r= (1, +m, c? + m, 2) (ec —c’ P,) P, |, 


= , a — , (fe — , i = 1 
Ba AChE) Pies ae (Canc P.) Ps |e P) Pe + (¢—e' P,) P,a (I, + mc?) (15) 


v 


+ my ¢, 1; (¢ — ec! Ps) P2l, + G 12) Be ah 


t = (¢—c' P,) P,a(e —e' P,) P, 1, | 
haben. 
Die in (14) angegebenen Werte der Eigenkreisfrequenzen sollten bei zweiachsigen Anhangern 


der behandelten Art beachtet werden, um im Bereich der zu erwartenden Fahrzustinde die Quer- 
schwingungsamplituden in zulassigen Schranken zu halten. 


4. Zusammenfassung. Fiir die Querschwingungen eines gelenkten Anhangers werden die beiden 
ungedampften Eigenschwingungszahlen und ein Verfahren zur Bestimmung des Einflusses der 
Dampfung angegeben. Bei geradliniger Bewegung des Zughakens wird ein vereinfachtes Anhianger- 
modell betrachtet. Die von L. Huber auf Grund von Versuchen angegebenen Ausdriicke fiir die 
auf einen Luftreifen wirkenden Seitenkrafte werden verwendet. Zur Bestimmung der beiden 
ungedampften Eigenschwingungszahlen werden kleine Querschwingungsamplituden betrachtet. 


(Eingegangen am 11. Januar 1958.) 


Anschrift der Verfasser: P. R. Paslay, General Electric Company, Schenectady, N. Y., USA. 
A. Slibar, Professor an der Techn. Hochschule Stuttgart, Stuttgart, Keplerstr. 10. 


1Siehe FuBnote 4 von Seite 384. 
? D. Williams, Proc. Inst. Mech. Eng., Autom. Div., 1951—52, S. 175. 
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Kritische Drehzahlen und Biegeschwingungen 
kontinuierlich besetzter Wellen bei verschiedenen Lagerungen* 


Von H. Bufler und H. G. Hahn 


1. Einleitung. Kritische Drehzahlen von Rotoren sind wegen ihrer technischen Bedeutung 
der Gegenstand zahlreicher Arbeiten!. Es liegt eine Reihe von Lésungen fiir verschiedene Be- 
setzungen der Wellen vor, wobei im allgemeinen eine zweifache kugelige Lagerung zugrundegelegt 
wird. 

Nachstehende Untersuchung bezieht sich auf den Spezialfall einer gleichmaig kontinuierlich 
besetzten Welle konstanter Biegesteifigkeit, fiir die im Fall der Rotation die kritischen Drehzahlen 
des synchronen Gleich- und Gegenlaufs, im Fall der reinen Biegeschwingungen die Eigenfrequenzen 
berechnet werden (denn in beiden Fallen liegt dieselbe Differentialgleichung vor) und zwar fiir 
verschiedene Randbedingungen. Mit dem Ziel, den Einflu8 der Tragheitsmomente auf die kri- 
tischen Drehzahlen und Biegeschwin- 
gungseigenfrequenzen bei den einzelnen 
Lagerungsarten beurteilen zu kénnen, 
wurde schon in einem fritheren Aufsatz 
der Verfasser? dasselbe Problem be- 7 
handelt, wobei man sich allerdings auf 
die jeweilige Grundfrequenz beschrankte. 

Die hier vorliegende Arbeit beinhaltet 
die vollstandigen Ergebnisse bis ein- 
schlieBlich der dritten Oberfrequenz und 
gibt Naherungsformeln fiir die weiteren 
(héheren) Frequenzen an, so daB das 7 
gesamte Frequenzspektrum erfaBt wird. NN A 

Folgende Lagerungsarten werden dis- 
kutiert (Abb. 1): 

a) FallI: Ein Wellenende .,einge- 
spannt*, das andere frei (flie- 
gende Lagerung), 

b) Fall IL: Beide Wellenenden kugelig 
gelagert, 

c) Fall 111: Ein Wellenende einge- 
spannt, das andere kugelig ge- 
lagert, 

d) FalliV: Beide Wellenenden ein- 


gespannt. 


Unter ,,eingespannt“‘ versteht man 
hierbei, da} die Welle zwar die Méglich- 
keit hat, zu rotieren, sich aber im be- 
treffenden Lager nicht senkrecht zu ihrer Abb. 1. Die untersuchten Lagerungsarten I bis IV. 
Achse verdrehen kann. 


2. Rechnung und Diskussion. Die hier benétigten Formeln kénnen z.'T. direkt aus unserer 
friiheren Arbeit? entnommen werden. Die eingeklammerten Zahlen vor den Formeln weisen 
auf die dort verwendete Numerierung hin. 

Es bedeuten I die Lange, E den Elastizitatsmodul, J das axiale Flachentragheitsmoment der 


Welle, m die Masse je Langeneinheit der besetzten Welle, @ das axiale Massentragheitsmoment 


* Vorliegende Arbeit entstand auf Anregung von Prof. Dr. H. Neuber, Lehrstuhl fir Technische Me- 
chanik der Technischen Hochschule Miinchen. ; 

1 K. Karas, Die kritischen Drehzahlen wichtiger Rotorformen, Wien 1935. 

2 H. Bufler und H.G. Hahn, Forsch. a. d. Geb. d. Ing. Wes., 23 (1957), 5S. 22. 
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je Langeneinheit der besetzten Welle, w die Kreisfrequenz bei Beriicksichtigung der Tragheits- 
momente, @, die Kreisfrequenz bei Vernachlassigung der Tragheitsmomente, ferner 


—@ bei synchronem Gleichlauf, 
O* =} 30 bei synchronem Gegenlauf, 
O bei Biegeschwingung, 


endlich «1 den Eigenwert bei Beriicksichtigung der Tragheitsmomente und («/)) den Eigenwert 
bei Vernachlassigung der Tragheitsmomente. 
Die Einfiihrung folgender dimensionsloser Gréfen hat sich als zweckmaBig erwiesen: 


; @* i 
(9) V= Pim?’ (1) 
a im b 55 je 2) 
(10) OTS EJ Zw. Wo = Wo EJ Py 
[2 2 O* 
(11) (Mae et (3) 


Es besteht folgender Zusammenhang zwischen den einzelnen Gréfen: 


(8) w@ [1 + V(al)?] = («I)*, (4) 
woraus, da fiir V = 0, w = w, und «1 = (a I)y ist, 

Dy = (0.1) (5) 
folgt. 


Das bezogene Tragheitsmoment V und das Frequenzverhaltnis w/w, lassen sich nun durch 
ol und p® ausdriicken. Es folgt namlich aus (4) unter Beriicksichtigung von (3), (2), (1) 


- ie 
=~ Gn f@le + py (6) 
Ebenso erhalt man aus (3) mit (6), (5), (2), (1) 
a — = = ae V(x 1)? =e joe . (7) 


Die fiir die Falle I bis IV geltenden Randbedingungen fiihren zu folgenden Beziehungen zwi- 
schen «1 und p*, welche die Eigenwerte reprasentieren (der untere Index bezieht sich hierbei auf 
die Lagerungsart): 


2 fg yes ae eee 
(16 a) s Nox l, p?) = = Venere sin (a l) Gin /(« 12 + p? 
pt + 2 (1? [(«l)? + p*] 


aneant py oC) GA Vot p+ 2=0, ®) 

(16b) OL p®?) = [2(« I)? + p?] sin (a 1) Gin /(@l? + p? = 0, (9) 
ro) _te(l IWVehr+p? _ 

Ce ne) ee ea ee oe 
(16d) Wo L, p?) = p® sin (wl) Sin Vw 1)? + p? 

+ 2a 1/ (012 + p? [1 — cos (x 1) Gof Va? + p?] = 0. (11) 


Die Gleichungen (8) bis (11) stellen je eine Schar mit unendlich vielen Kurven dar. Zu jeder 
kritischen Drehzahl bzw. Eigenfrequenz als Funktion der Tragheitsmomente gehért eine be- 
stimmte Kurve', was durch den oberen Index y angedeutet sei (also z. B. y = 1 fiir die Grund- 
frequenz, wofiir the == 0, ie = 0... zustandig sind). Diese Kurven wurden numerisch fir die 
Grundfrequenz (vy = 1) bis einschlieBlich der dritten Oberfrequenz (vy = 4) ausgewertet und in 
Abb. 2 dargestellt, wobei sich die verwendete logarithmische Auftragung als zweckmabig erwies. 


* Streng genommen sind es je zwei Kurven, die zur p?-Achse spiegelbildlich liegen und beide zum selben 
Ergebnis fiihren, 


y 
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Zu ihrer Diskussion seien einige Bemerkungen gemacht. Fiir p? = co haben samtliche Kurven 
zur p*-Achse parallele Geraden als Asymptoten, fiir p? <0 sind die Kurven nur im Bereich 


—(«l?<p?<0 reell. In den 
Tabellen | bis 3findet man die « I- 
Werte fiir p? = 0 sowie die Asym- 
ptoten fiir p? = co und die Punkte 
p? =—(«l)? nebst den dazuge- 
hérigen Bestimmungsgleichungen. 
Man kann daraus entnehmen, da 
im Bereich héherer Frequenzen 
(v groB) die maximale Abweichung 
Axl der einzelnen «1 bei zweifa- 
cher Lagerung der Welle (Falle II 
bis IV) einen konstanten Wert 
annimmt, der relative Unter- 
schied Ax I/x 1 mit wachsendem y 
immer mehr verschwindet. Die 


ron fs y oe 
Kurven f,, = 0 und pan = 0 kén- 


nen daher fiir groBe » praktisch 
durch die Kurven J i so Mids hi. 
sin (¢ 1) = 0 baw. al=va (vy=1, 
2,...) ersetzt werden, welche 
zur p*-Achse parallele Geraden 
sind. Ebenso laBt sich zeigen, daB 


fiir groBe y die Kurven f\” = 0, die 
fiir die einseitig eingespannte Welle 


gelten, durch a! =~ (2y—1) 
ersetzbar sind. - 

Um das Frequenzverhaltnis 
@/@,) als Funktion des bezogenen 
Tragheitsmomentes V darzustel- 
len, braucht man nur die in 
der «al-p?-Ebene  dargestellten, 
die Eigenwerte  beinhaltenden 


Kurven fas 10. ie 0; iat ==) 


und i = 0 mittels der Trans- 
und (7) 
in die Zz V) Ebene abzubilden. 
Dabei ein aseten sich die in 
Tabelle 4 aufgefiihrten Punkte. 
Man entnimmt aus der letzten 
Zeile der Tabelle die dem Gleich- 
lauf gesetzten Grenzen. Da hier 
V =—O/lm ist, folgt, daB 
im Fall I fiir 

) 4 

Pia ~ m2 (v— 1p” 
in den Fallen II, III, IV fir 

(0) 1 


im 


formationsformeln (6) 


y? 7 


Tabelle 1. Spezielle Werte der Kurven ft) = (). 


Bedingungen pi 0 | 


Bedingungs- 
gleichungen 


y = 
2 4,094 
3 7,855 | 
4 10,996 


cos (al) = 


lean poloo wlQ 
a 


to| 2p) 


p?> 0 p? = p* = —(«1)? 


| cos(al) = 0 


wla wlan wile wa 
i 


Bedingungen 


Pp SO 


Bedingungs- 
gleichungen 


Pee ll coll == SS PAl 
2 7,069 
3 10,210 
4 13,392 
A | 
n = (" + ral | 


m7 


tg(al)—al = 0 


al = 4,493 
7,724 
10,905 
14,064 


Bedingungen p?=90 poo p? = p™* = —(al)? 
3edingungs- Kee ea eae 
Bedingwnes- | gof(at) cos(at)—1—=0| tg(a=o | oe tay 2-0 
ell l= 4.730 | al=x al = 6,284 
2 7,853 | 2% 8,991 
3 10,996 37% 12,568 
4 14,137 A xt 15,453 
| . 
n ~(+5 Tt na =(n-+ 1)a 


kein synchroner Gleichlauf méglich ist. Bei einer zweifach gelagerten Welle ist die Grenze dafir 
somit unabhangig von der Art der Lagerung. Auf diese Tatsache hat bereits Grammel' hingewiesen. 


1G. B. Biezeno und R.Grammel, Technische Dynamik, 2. Bd. S$. 201 u. 202, Berlin/Géttingen/Heidel- 


berg 1953. 
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Tabelle 4. Zuordnung einiger Punkte der — - V-Ebene zur « |-p*-Ebene. 
0 


| 
Pp? | (a 1)? [4 w/o 
i 
De; met oe | 0 Bereich des Gegenlaufs 
0 (« 1)? 0 1 und der Biegeschwingungen 
D0 
1 : 
ym bei II, II,1V = bel dT dy 
ba ee es Bereich des Gleichlaufs 
: l , 1 ; 
( — A cre Deianal =e Se - bei 


Die(—s V\-Kurven wurden numerisch ausgewertet und in Abb. 3 bis 6 fur die Falle I bis IV 
Mo 


graphisch dargestellt. Da bei den héheren Frequenzen (y groB) die Kurven fo ==10, fee ==) 


naherungsweise durch f= sin (xl) =0 (daraus «1 =yz) ersetzt werden diirfen, verschwindet 
hier der bei den tieferen Frequenzen nicht unwesentliche Einflu8 der Lagerungsart und zwar um so 
stirker, je gréGBer y ist. Man kann dafiir die fir Fall II streng giiltige Formel anwenden, die man 


aus (6) und (7) mit al =yz erhalt: 


Dare i 1 . (12) 
% YV1l+wapPpV 
Man sieht in der Tat aus Abb. 4, 5 und 6, daB fiir vy = 4 nur noch geringe Unterschiede zwischen 
den Fallen IT, III und IV bestehen, so daB man fiir y > 4 ohne weiteres von der Naherungsformel (12) 
Gebrauch machen kann (die in Abb. 5 und 6 vergleichsweise gestrichelt eingezeichneten Kurven 


wurden fiir » = 4 nach (12) berechnet). 
vA 
| | 
b | y 
U 


Abb.7. Die im Beispiel behandelten Wellenbesetzungen 
a) R/L=0,1,  -b) R/l = 0,3. 


~+—2R—e| 


2R 


Ebenso darf man im Lagerungsfall I fiir groBe » naherungsweise 
f” = cos (a1) =0; daraus = a 1 = = (2 —1) 


setzen und bekommt wegen (6) und (7) 
a) 1 


= —— ee a . (13) 
i= 2 
/ He B (Q»— | V 

Die danach fiir y = 4 berechnete Kurve ist in Abb. 3 gestrichelt eingetragen und zeigt eine etwas 


gréBere Differenz mit der genauen Kurve als bei den anderen Lagerungsarten. Da diese jedoch mit 


wachsendem y schnell abnimmt, stellt (13) fiir » > 4 eine sehr brauchbare Naherung dar. 
Damit ist praktisch der gesamte Frequenzbereich erfaBt. 
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In Abb. 3 bis 6 sind die fiir den Gleichlauf, den Se balla 
Gegenlauf und die Biegeschwingung geltenden ee tpl 
. . ER pee eae en 
Relationen zwischen V und @/I?m angegeben. a ee 
AuBerdem sind die Formeln fiir die Frequen- | oa 
h Beriicksi hti cs ; | | Stir on | 
zen Wy ohne Beriicksichtigung der Tragheits- es Sachin etalon 
momente in den einzelnen Diagrammen an- | Pe Pen eee 
gefiihrt, so daB man bei gegebenen Daten e ea ees | 
sofort das Ergebnis fiir die kritischen Dreh- ek zeae 
ee ee 
zahlen ie aga); CF 
30 /w a lye ee eh 
a pe (me KE (in U/min) eo mao) oo 
0 3 | 
5 . cae 
oe . . . . sv 
baw. fiir die Biegeschwingungseigenfrequenzen Pe. ee 
a4 oe | AS 
ee : Ces 16 60 | ot 
ES Sesilag on AE) fi] | See |e | 
0 he 
ale =a 
angeben kann.. Die Formel fiir «, folgt aus (2) | ead eee ta ae Oe eee 
| * ateet 
und (5) und lautet ee i a ea racer ear er ae 
——— S ee titer iealleee cutitaa lob eke a 
anh) te. | eu (eee Re ec 
Os = (x 1); an . (14) = | ~| ARS faa 
as | | i] ob. oe OD De 
3. Beispiel. Es ist eine masselose Welle der 2 eee | 
- . . . . . 2S 
Lange /, die mit dicht aufeinandersitzenden © | | See | SEE 
. ° i 
Scheiben vom Halbmesser R besetzt ist, ge- 3 <| | es 
3 a 2 Ry Sq9e |aqa 
geben. Das axiale Massentragheitsmoment ‘3 | 3 eee (eee 
der Scheiben sowie ihre Masse je Langeneinheit = 
) 
. aD = ne a a 
seien bekannt (0 und m). Es soll fiir den 3 [6 ls 42 Loge Pea ate | On ga ape 
synchronen Gleich- und Gegenlauf und firreine S| 2 8 _ | J aa l|as|/|as) | as 
: : roe Blea ales ae cS to je 
Biegeschwingungen das Frequenzverhaltnis & : eas ook Hb Sos tier 
©/@, ermittelt werden und zwar fiir die Lage- $ BESS) or | ee ene ae Cen meer 
e . . ae Lo. | ; 
rungsfalle I bis IV, wobei das Verhaltnis R/l > = J) ASS ASS SSS (RSS 
zu 0,1 bzw. 0,3 (Abb. 7) angenommen sei. Es | Saou eee eS aaa et 
gilt dann z e 
6 1/R\ 0,0025 , a e |i 
Sat SS | |S) Eta 
Pim Amey | baw. 0,0225. 8 £) a 
. : ; 5 Sea] 
Mit Hilfe von Abb.3 bis 6 erhalt man 8 aS 
. . ~ ima 5 
sofort die in Tabelle 5 zusammengestellten S (eae 
Ergebnisse. = 
fia HC : eye 19 ge is 
4. Zusammenfassung. Fiir eine gleichmaBig 2 | 3 
kontinuierlich besetzte Welle konstanter Biege- o7 
“ce . . elie is} St 
steifigkeit werden die kritischen Drehzahlen = ees 
des synchronen Gleich- und Gegenlaufs und die ere | 
Eigenfrequenzen der Biegeschwingungen be- | oe 
rechnet, und zwar fiir vier verschiedene Lage- 5 Ee eet eee es 
rungsarten (fliegend gelagerte Welle; zweifach q psillies a) | Mache ee |e Rte |e Ses 
. A * fale ets exe ela Senet One a Spee 
kugelig gelagert; ein Wellenende,,eingespannt“, eT a Ia ae get hace acca Rigs 
aq s . | eed = ; = 
das andere kugelig gelagert; zweitach ein Be Al eee Pee cpianee tn heres (cea 
gespannt“‘). Die genauen Ergebnisse fiir die '& [ll Shs SERISRRLSRS 
Grundfrequenz bis einschlieBlich der dritten | espace RT ee eee ae ee 
Oberfrequenz sind graphisch dargestellt und : ; ; : 
: 2 ; a copay eS ES 
lassen den EinfluB der Tragheitsmomente und 2 3 ts ee Bare E A| as 5 
at SSO SS oe eres oS 
der Lagerungsart erkennen. Fiir die héheren ad 2 ad 2 a eeia'a8 
es . a ne i ey a 
Frequenzen werden Naherungsformeln mit- ‘BOS |S OP) Ss oP |S we 
Oom)/oOoOmM;/o0Mm |;C°oOn 


geteilt, so daB das gesamte Frequenzspektrum 
erfaBt wird. 


an a a a 
(Eingegangen am 16, Januar 1958) P = e ad 2 
B= ae | 3c — 
Anschrift der Verfasser: Dr.-Ing. Hans Bufler und ss = 54 - SS = os a 
Dipl.-Phys. Hans Georg Hahn, Miinchen, Technische 8 3 8 3 


Hochschule, Institut fiir Technische Mechanik. 


** Es gelten jeweils dieselben Ausdriicke wie im Lagerungsfall IT. 


* Diese Werte liegen auBerhalb der Abb. 5 u. 6, kénnen aber leicht mit Hilfe der Abb. 2 ermittelt werden. 


*** Fir »y >7 tritt hier kein Gleichlauf mehr auf. 
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Berechnung von Biegeschwingungen abgesetzter Wellen 
mit Zwischenbedingungen mittels Ubertragungsmatrizen 


Herrn Professor Dr. A. Walther, Darmstadt, zum 60. Geburtstag 
Von R. Zurmihl 


1. Einleitung. Kompliziertere Aufgaben der Biegeschwingung, insbesondere Schwingung ab- 
gesetzter Wellen mit Zwischenbedingungen, aber auch ebener und raumlicher Stabwerke sind in 
neuerer Zeit mehrfach mit Hilfe von Matrizen behandelt worden!. Dabei wurde die Rechnung, aus- 
gehend von bekannten Randbedingungen am Balkenanfang, mit probeweise angenommenen Fre- 
quenzwerten durchgefiihrt, die solange abgedndert werden, bis die am Balkenende geforderten 
restlichen Randbedingungen erfillt werden (RestgréBenverfahren). Dieses Vorgehen eignet sich 
besonders fiir den Einsatz digitaler Rechenautomaten, da der gleiche und in sich sehr homogene 
RechenprozeB viele Male wiederholt wird, was der automatischen Rechnung sehr entgegenkommt. 
Fiir die Handrechnung ist indessen eine unmittelbare Frequenzberechnung wiinschenswert. Im 
folgenden wird ein gleichfalls auf der Verwendung sogenannter Ubertragungsmatrizen beruhendes 
Naherungsverfahren entwickelt, das die Aufgabe der Biegeschwingung abgesetzter Wellen stiick- 
weise konstanter Querschnitte mit beliebigen Rand- und Zwischenbedingungen in eine gewohnliche 
Matrizeneigenwertaufgabe iiberfiihrt, die sich durch gute Konvergenzeigenschaften zur iterativen 
Bestimmung der niedrigsten Eigenfrequenzen des Systems auszeichnet. AuBer den Kigenfrequenzen 
erhalt man die zugehérigen Schwingungsformen ohne zusatzlichen Arbeitsaufwand. 


2. Herleitung der Matrizengleichungen. Die auf Biegeschwingung zu untersuchende Welle 
(Balken) sei in n Abschnitte (Felder) von jeweils konstanten, nicht notwendig verschiedenen Quer- 
schnitten unterteilt. Anfang und Ende des i-ten Ab- 
0 Bre © He @ schnittes der Lange /;, der Biegesteifigkeit E J, und 
a ee der Massenbelegung yw; je Langeneinheit seien mit 
0 1 ede HEY Se en ss) n t—1 bzw. i bezeichnet (Abb. 1). Feldgrenzen sind 
Abbe ls Bezeichnune der Balkenabscharttc (starre oder elastische) Stiitzen, Gelenke, Hinzel- 
massen, Querschnittsspriinge, aber auch Stellen 
bloBer Unterteilungen, die aus Griinden der Genauigkeitssteigerung vorgenommen werden. Unter 
TU Rpanoae ae der Schubverformung gilt fiir das i-te Feld die Differentialgleichung der Biege- 
inie 


E 1, w(x) = q(x) (1) 
mit der Belastungsdichte q(x), fiir die im Falle freier Biegeschwingung die Tragheitslast 
q(x) = (4, 0 w(x) (2) 


mit der — einstweilen noch unbekannten — Eigenfrequenz @ einzusetzen ist. Fassen wir die Be- 
lastung q(x) zunachst als gegeben auf, so ergibt sich als Lésung von (1) ein linearer Zusammenhang 
zwischen den Zustandsgréen — Durchbiegung w, Neigung w’ (bei fehlender Schubdeformation 
gleich der Querschnittsneigung —q), Biegemoment M = —E I, w'’ und Querkraft Q = M’ — in 
Feldanfang und Ende in der Form 


2 I: 


HO Wie alg Wee by edly Q._1¢ Ww; , 
/ I 1; I? 
a wi-1 —M;_1—=— —Q;-1>-5-+ Hee| 
E I; ed) de v ( (3) 
M;= Vie SG ers + M*, 


C= OF es sls 


1 E. Pestel, Abh. Braunschw. Wiss. Ges. Bd. 6 (1954), S. 227; S. Falk, Abh. Braunschw. Wiss. Ges. Bd. 7 


(1955), S. 74; H. Fuhrke, Ing.-Arch. 23 (1955), S. 32 » : 
PORE A ey ARES Oe oe ), S. 329, Ing.-Arch. 24 (1956), S. 27; H. Unger, Intern. Koll. 
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mit den Belastungsgliedern! 


O(a) = — fale ae, oe 
M*(x) =f or) dé , M* = M*(i) , | 
eee (4) 
w(x) == EE | M*(€) dé, Wer ab *(1,) | 
w*(x) =f w'*( dé, Fico ALIN 


Die Vorzeichen sind dabei nach Abb. 2 gewahlt. 
Fassen wir die vier Zustandsgré8en zum Zustandsvektor 3,, die entsprechenden Belastungs- 
gréBen zum Belastungsvektor 3; zusammen: 


Q 
ia et w;* (5) u{ ——— by 
Pu @ 
L 


(i M,|’ ee cal 
(o,) (os) ole 


so schreibt sich (3) in Matrizenform Abb. 2. 
Vorzeichen der Zustandsgréfen. 
3: = Bi 81-1 + 37 (6) 


mit der von den Feldabmessungen abhangigen Feldmatrix = Ubertragungsmatrix 


rp 3 
(lt —s5 —sz7\ 
I RP 
Er 2EI\: (7) 
0 1 l 


0 
Voom to Loge 


Die Feldmatrix 7%, tibertragt die am Feldanfang herrschenden, im Vektor 3, _, zusammengefaSten 
ZustandsgréBen auf die neuen, am Feldende herrschenden unter Hinzutreten der Lasteinfliisse 3. 


%,=|9 1 


3. Naherungsformel fiir den Belastungsvektor. Die Komponenten des Belastungsvektors 3* 
errechnen sich nach (4) aus der Feldbelastung q(x), die zufolge (2) der unbekannten Durchsenkung 
w(x) proportional ist. Diese aber ]aBt sich angenahert mit Hilfe der Werte w; an den Feldgrenzen 
ausdriicken, und zwar bedienen wir uns dazu zweckmaBig einer sogenannten Hermiteschen Inter- 
polation, bei der als Stiitzwerte auBer Funktionswerten w; auch noch Ableitungen bis zu einer 
bestimmten Ordnung benutzt werden, die ja im Zustandsvektor bis zur dritten Ordnung vorkom- 
men. Einfachheitshalber beschranken wir uns hier auf Mitnahme der ersten Ableitungen und be- 
nutzen als Stiitzwerte der Interpolationsformel fiir w(x) des i-ten Feldes die vier Randwerte w;_1, 
W;—1, W;, W;, nahern also w(x) durch eine — recht genaue — kubische Interpolationsparabel an. 


Bezeichnen wir voriibergehend die Feldgrenzen mit 0 und 1 und setzen die Feldlinge J; gleich 1, 
so lautet der Naherungsausdruck fiir w(x): 


w(x) = Lo(x) wy + Ly(x) wy + Ho(x) wy + H,(x) w, , (8) 
wo sich die vier Interpolationspolynome dritten Grades I,(x), H,(x) aus den Forderungen 
£0) =1, - 1,0) =0,- Li(0)=0, LG) =90, 
OMS URL (Leal, Poe (0) 3052" (1) 0; 


1 Nach S. Fall, Ing.-Arch. 24 (1956) S. 216. 
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ergeben zu 

L(x) =1—3 4 4 2 4; 

L(x) = 3 2—2 2%; (9) 
H,(x) = x — 2 «7 + 23, 

HA, (x) = —x#+ 23. 


Hiermit errechnen sich dann die Belastungsglieder aus (4), wenn wir noch wieder auf die Feldlange 


I], = 1 tbergehen, zu 


JI I f / 
OF = 0815 Geir bm) yy (ala ay] — m0 a, 
l l f / ‘ 
M* =—po' P E (7w;—1 +3w,) + 60 (3 w;—_1 —2 rm) + 0,0 w;,, a5 
Pe Rae haa! zh TS ee 
Dips ae ky (4 w;—1 + w;) + 130 (2 wj3-1— w,)I, 
Fy ef 


Ww; = Sa pie 


if 1 / 47 
Fe (6 Wi—1 ++ w;) -|- 1260 (5 Ua Co x w)| . 
Dabei ist ttherall ] =1,, EI = EI,, u =; 20 setzen. Zugleich haben wir hier noch den EinfluB 
etwaiger diskreter Einzelmassen m; und Einzeldrehmassen @, in Q; und M, beriicksichtigt, wobei 
diese Massen jeweils am Feldende sitzen. 


4, Einfiihren bezogener Grofen. Zur Vereinfachung der Zahlenrechnung fiihren wir bezogene 
ZustandsgréBen ein unter Verwendung einer — an sich beliebigen — Bezugslange /, einer Bezugs- 
steifigkeit E I und fir das weitere einer Bezugsmassenbelegung uu. Als endgiltige ZustandsgréBen 
wahlen wir, unter Vorzeichenumkehr bei w zur Vermeidung negativer Elemente in %, 


w=—w, 9 =—w'l, M=—, Q = <— (11) 


— 


als Komponenten des Zustandsvektors 4,. Mit den Verhaltniszahlen 


li El; 


B. cay) i 
Be a (12) 
wird die Feldmatrix 
3p2 ps 
| 
68 uo 38? ii 
Sosa | nese ls (13) 
@ I pb 
lo 0 0 ] i 
Mit den weiteren Verhaltniszahlen 
Lai eahe ae twa lee Mr iae (as) 
sowie den Frequenzparametern 
cons iw lt es ja ‘ 
yi ee A= 3550 (15) 


erhalten wir dann fiir die Komponenten des Belastungsvektors 3; aus (10): 
= 4 _ a Fy ™ 
ws = Ase (90 w;— 1 + 106 gi_1 els Wa AB) » 
—, 40/3 = 7% 1 p 
pi = A~— (336 wi_1 +428 G;:1 + 840,—218G)), (10a) 
M3 = A 0f?(147w;_1 +218 9;,_1+ 63 w, — 148 %,) — 420 Az; ¢;,, 


QF = A of (210 w;_1 + 358 G,_1 + 210 w, — 358) + 420A.6, i, 
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wo iiberall vereinfachend 0, B, « anstatt o,, B,, «; geschrieben ist. Fassen wir die hier auftretenden 
Biege-Stiitzwerte w,, ~, zum Teil-Biegevektor 


W;—1 
ee | le (16) 
Ww. 


\5, 


die Koeffizienten in (10a) aber zur Belastungsmatrix 


900 f/x l0ef ix 150 Bla —4o f®/a 
% = 3360 A/a 420 f/x 840 f/x — 21 9 B4/x (7) 
‘1147062 210f® 63062  —14(0 83 +307) ; 
21006 = 350f? 210(08+2c) —35ef? }; 


zusammen, so erhalten wir als Naherungsausdruck fiir den Belastungsvektor 
4, = 1B; 0; (18) 


Damit lat sich nun die Rechnung, zunichst fiir den Balken ohne Zwischenbedingungen, auf 
folgende Weise durchfiihren. 


5. Der Rechnungsgang. Man beginnt die Rechnung am Balkenanfang i = 0, wo von den vier 
ZustandsgréRen des Vektors 4), den wir in der Form 


{I 0 0 f 

a 1h) TEE O.|| 5 One 

60 =} 9 ore 0 Po + 1 Moa 0 Qo (19) 
0 0 0 1 


aufspalten kénnen, zwei aus Randbedingungen bekannt sind; z. B. ist bei freier Lagerung 
Wy = M, = 0, womit sich (19) auf 


0 

0 —a 

; Qo (19a) 
1 


Se 
—) 
| 
coor oS 
Ss! 
i=) 
+ 


reduziert mit zundchst unbekannten Faktoren qo, Q,, welche durch die gesamte Rechnung als 
Unbekannte mitzufiihren sind und sich zum SchluB aus den am Balkenende gegebenen zwei wei- 
teren Randbedingungen ermitteln lassen. Allgemein wird der Anfangsvektor 3) die Form 


jo = Cy Xy + Cy Xe (20) 


haben mit zwei Unbekannten X,, X, und zwei Vektoren ¢,, ¢,, die im einfachsten Falle, wie oben, 
Einheitsvektoren sind, davon aber auch, wie sich noch zeigen wird, abweichen kénnen. Dem Vek- 
tor 3) sind somit zwei Spalten zugeordnet, die man sich mit den beiden Unbekannten X,, X, als 
Faktoren versehen denkt, und die nun der fortgesetzten Multiplikation mit den Feldmatrizen ay 
der aufeinander folgenden Balkenabschnitte unterworfen werden. Dabei bedient man sich zweck- 
maBig einer von Falk vorgeschlagenen Anordnung fiir die Multiplikation von Matrix mal Vektor 
nach Rechenschema 1. — Zu diesen beiden X-Spalten treten nun der Reihe nach weitere, von den 
Belastungsvektoren 3* herriihrende Spalten, versehen mit den gleichfalls noch unbekannten Fak- 
toren A w,, A @;, und zwar fiir jedes Feld ein Paar neuer BiegegréBen, sofern nicht, durch Rand- 
und Zwischenbedingungen, einige dieser GréBen zu Null werden, wobei die zugehérige Spalte ent- 
fallt. Die Vektoren 3, erscheinen somit in eine wachsende Anzahl von Spalten zerlegt. Auf alle 
diese Spalten wird fortgesetzt die Multiplikation mit der F eldmatrix, 7; 3:1 angewandt, wobei 
jeweils in den vier letzten Spalten noch der nach (18) gebildete Lastvektor 3* als Summand hinzu- 
geschlagen wird. Auf diese Weise baut sich ein Rechenschema 2 auf, in dem auBer den als Faktoren 


97% 


ad 
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fungierenden Feldmatrizen %; auch die Belastungsmatrizen 8. eingetragen sind. 


Rechenschema 1 


Am Balkenende, i = n, treten nun zwei weitere Randbedingungen auf, z. B. bei freier End- 


lagerung w, = 0, M,, = 0. Dies ergibt zwei Gleichungen fiir die beiden Anfangs-Unbekannten X,, 
X,. Fassen wir diese zum Vektor r, die BiegegréBen w,, 7, aber zum Gesamt-Biegevektor ty 2u- 
sammen: 


ii) =| |I (21) 


so schreiben sich die Endbedingungen in der Form 


Kr+Adw=0, (22) 


wo die Koeffizienten der Matrizen ‘8, 9, unmittelbar den beiden Zeilen von 3,, welche den End- 
bedingungen entsprechen, zu entnehmen sind. Ist ‘3 nichtsingular, was die Regel sein wird, so 
lassen sich hieraus die beiden Unbekannten ¢ in der Form 


(23) 


C= ip 
Rechenschema 2 zum Aufbau der Zustandsvektoren 3, 
| zt Aso | 
op oo se re = = = | aa 
2s XxX, Aw, | AM Aw, | AQ, | Aw, | AGs | cag An | Agyn | 
! a ee ea ee 
| | oa Geer | 
| | 
| 30 
| | 
Soieta| 3a ee ae 2 |_| 3 
at lod 
¥, O1 | | | 
| 2 
| 01 
| | | 
Pe ee Lat “ Siceen| 2 ie | | 
| | ie | |" 
B, Se | | 
| | | 
| | | | be 
aise a es er oa pet Z ate SS, Ml | 
ay msi Leds | st | ll | re | EP 
V C | 
Bn Wn | 
| | 3 
| bn 


durch Elimination ermitteln, wobei die Matrix © sich als Lésung des Gleichungssystems 


FE+Q=0 
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ergibt, wie leicht nachzurechnen. Hat man auf diese Weise die Unbekannten y in Abhangigkeit 
von v ausgedriickt, so ergeben sich samtliche BiegegréBen tv aus den jeweils ersten beiden Zeilen 
der Vektoren 3;. Sie erscheinen im Schema in der Form 


w=Rr+ASGw, (25) 
woraus durch Einsetzen von (23) fiir die BiegegréBen folgt 
pw =ARC + S) vw =AAw (26) 
oder mit 
1 2920 
ao eae 
endgiiltig 
Win = x w (27) 
mit der Matrix 
AW=RKC+ 6S, (28) 


die sich unmittelbar an Hand des Rechenschemas 3 aufbaut. Damit ist die Aufgabe der Biege- 
schwingung auf eine Matrizeneigenwertaufgabe zuriickgefiihrt. 

Die numerische Lésung der Eigenwertaufgabe (27) geschieht am einfachsten iterativ nach dem 
bekannten v. Misesschen Iterationsverfahren!, das bekanntlich gegen den hetragsmaBig gréBten 
Eigenwert x, konvergiert. Dem entspricht der kleinste Wert A,, der @? proportional ist, dem 
Quadrat der in der Regel interessierenden niedrigsten Eigenfrequenz. Da bei Biegeschwingungen 
die Eigenfrequenzen @; mit steigender Ordnung j angendhert wie die Quadratzahlen j? zunehmen, 
die Werte A, also wie j*, so ist im allgemeinen sehr gute Konvergenz des Iterationsverfahrens zu 
erwarten, was durch Zahlenrechnungen bestatigt wird. Die beim Iterationsverfahren gleichfalls 
anfallenden Eigenvektoren {v stellen zugleich die Biegeschwingungsformen in Gestalt diskreter 
Ordinaten nebst Steigungen dar. Hiernach la8t sich dann auch der Ausgangsvektor des Iterations- 
verfahrens hinreichend gut schatzen, so daB sich die Rechnung rasch einspielt. 


Rechenschema 3 
sur Aufstellung der Matrix YA 


xt | Ap | 
~ | 
| @& | (&) 
| | 
| ®& 
S | Elimination | 
Ss 
etl 
1 0 
(v’ 
0 1 


6. Beispiel. Fiir die abgesetzte, beiderseits frei gelagerte Welle nach Abb. 3 mit 


l=T, EI=EL, ii =, 4ET ET 


in an 
o = 4, 0=2, py =1, 


X& = 1, Onl (ia ie ii et sell 


Abb. 3. Abgesetzte Welle 


1 Vgl. z. B. R. Zurmihl, Praktische Mathematik fiir Ingenieure und Physiker. S.159ff. 2. Aufl. Berlin, 
Gottingen, Heidelberg 1957. 
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erfolgt die Rechnung nach dem Schema 4 


Rechenschema 4 


—_—_—_——— ee ca ea 
Po Qo | AM Aw, AQ, A 2 
ee  ——————e  — — — — _———__eesw 
| 0 @ | w, = 0 
1 Ona Po 
0 0 
BS; Ci 0 1 
As aS 7B rae i sdess/ Ane A 1 1/4 5 ee) io, 
LNG a2 —10,5 0 1 “3/2 3/4 1 3/4 21 42 — 10,5 1 
294 42 126 —28 0 0 1 0 1 42 126 — 28 
420 70 420 —— 0 WOO) a0 iI 0 1 | 710 420 — 70 
3 90 10 15 ; —— A ih ais! 1 2 5 222 937,5 —156,5 — 4 0 
336 42 84 =e eee) 3 1 39/4 483 2394 —346,5 —21 G2 
L472 S63: —14 O QQ i il 0 2 ial 693 — 77 —14 0 
210 35 210 shy Wi) Ww 1 0 1 70 630 — 35 —35 
: ; 2 5... 5202» 087 sat —156s5 feet 
| 0 2 e2 693 — 77 —14 
| = —— | 
29 — 56 1 AG, 
a 397,5 —346,5 | 1 | Aw, 
18 38,5 1 AG, 
ahh ial ay 1 A Qs 
; 29. ,397,53' —18 15.5, o)lseees 
9 20 318,375 —10,375 —13,75 aD 
8 179,625 0,375 —10,25 | =9, 
+34°--586,875 10,873)” 3170 (ee 
Die Durchfiihrung der Iteration an der Matrix 8 2{ zeigt das Rechenschema 5 
Rechenschema 5 
Po wy MA Pe 
232 3180 —144 —124 
160 2547 — 83 —110 
64 1437 3 — 82 
—272 —4695 87 254 
2, 2 1 a 3 
7052 5661 Syst —10609 
1,25 1 0,57 —— Megexyy 
3619,8 2905,4 1672 ape 5 460,4 
1,2459 1 0,5755 i sane SMR TOR 
3619,22 2905.31 1672,57 — 5461,18 
1,245 726 1 0,575 694 — 1,879 724 
3619,194 2905,303 1672,591 — 5461,202 
1,245 720 1 0575703 | —— 1,879736—= wy, 
2905,303 | | 
8 x, = 2905,303 
pte 53630163 
i$ = 6,939 03 
i,= 1,623 02 
exakt 1,620 43 
Fehler 0,002 59 = 1,6%/o 


Zur Berechnung des nachsthéheren Eigenwertes A, kénnte man die Iteration mit zum ersten Eigen- 


vektor {v, orthogonalen Vektoren durchfiihren (Verfahren von Koch); doch ware dabei feinere 
Unterteilung der Welle anzuraten. 
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re Stiitzfedern. An durch Federn oder Drehfedern elastisch gestiitzten Balkenstellen treten 
Spriinge in den Querkraften bzw. Biegemomenten auf, die der Durchsenkung w, bzw. Neigung 9; 
proportional sind. Fiihren wir an Stelle der Federkonstanten c; bzw. Drehfederkonstanten C, di- 
mensionslose Zahlen = = 
ae Tie! 

6EI ao Fad 
ein und bezeichnen die Stelle unmittelbar vor der Stiitze mit i’, die unmittelbar dahinter wieder 
mit t, so ergeben sich die Sprungbeziehungen 


Vi (29) 


: Q; 3 Qi meat w; ? M, = My = L, Pi z (30) 
Der Ubergang vom Zustandsvektor jv vor auf den Vektor 4; hinter der Stiitze i 1aBt sich dann gleich- 
falls in Matrizenform schreiben zu 
8: = Bae (31) 
mit einer sich auf den Stiitzpunkt i beziehenden Ubertragungsmatrix, einer Punktmatrix 
. Ie Cr 60.10 
a tall Onnilin EOia0 
Eracltsd re eet @) 
—-y, 0 01 
An Stelle von (6) tritt dann 
i= BOs bai Be - (6a) 


Die hier auftretende Produktbildung $f, 3; von Punkt- und Leitmatrix lat sich nun nach einem 
Vorschlag von Falk auf vereinfachte Weise vornehmen, indem man an Stelle der Feldmatrix eine 
nichtquadratische, aus vier Zeilen und sechs Spalten bestehende sogenannte Leitmatrix &, ver- 
wendet, welche Feld- und Punktmatrix in sich vereinigt. Sie lautet hier 


sei ee 0a 0 
o o 
6B 3B, 
Gen mer ee pee oral 
OC -=0.5 0 be = On). 
Damit schreibt man dann anstelle von (6a) 


bi = Us bi-1 + 87> (6b) 
wobei aber nun abweichend vom Ublichen die Matrizenproduktbildung so verstanden sein soll, 
daB die erste bis vierte Spalte von &, wie iiblich skalar mit den vier Komponenten von 3;_,, die 
beiden letzten Spalten aber mit den beiden ersten Komponenten w, und @; des neuen Vektors 3; 
multipliziert werden. Die Produktbildung erfolgt also iberlappend. Das ist dadurch méglich, 
da8 die Zusatzspalten von &, in ihren beiden ersten Zeilen, welche die neuen Werte w; und @; erst 
ergeben, Nullen aufweisen. Die Zusatzspalten treten also erst in Aktion, nachdem diese beiden 
neuen Komponenten von 3; bereits vorliegen, nimlich bei Berechnung von M; und Q,, die dann 
die Werte rechts vom Stiitzpunkt, also mit Zuschlag der Spriinge darstellen. Die Werte links vom 
Stiitzpunkt ergeben sich ohne Beriicksichtigung der Zusatzspalten der Leitmatrix. — Bei Federn 
am Balkenanfang sind schon die Anfangsvektoren fiir die Unbekannten w, baw. Py entsprechend 
abzuandern in 


iy | Go 
i] 0 
0) 1 
0 | I a 
el 


8. Innere Randbedingungen (Zwischenbedingungen). AuBer den an den beiden Balkenenden 
gegebenen je zwei Randbedingungen kénnen auch an Zwischenstellen geometrische oder dyna- 
mische Bedingungen als innere Randbedingungen auftreten, was das Auftreten zusatzlicher Un- 
bekannter zur Folge hat. Die praktisch bedeutsamsten Faille sind: festes Lager mit der Bedingung 
w, = 0 und Gerbergelenk mit der Bedingung M; = 0. Zusatzliche Unbekannte ist im ersten Falle 
die Lagerkraft A, gleich einem Querkraftsprung AQ;, im zweiten Falle eine Winkelanderung Ag;. 
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Der Rechenablauf wird hiervon nur in der Weise beriihrt, als auBer den beiden Anfangsunbekannten 
X,, X, noch weitere Unbekannte X,, X,... auftreten, die in die Rechnung eintreten, sobald man 
an die Stelle der inneren Randbedingung angelangt ist. Es erhdht sich also lediglich die Zahl der 
Unbekannten-Spalten, die sich nach und nach auffiillen, und der Vektor y der Unbekannten X; ent- 
halt mehrals nur zwei Komponenten. Ist die Rechnung am Balkenende angekommen, so hat man auBer 
den beiden Endbedingungen nun auch noch alle Zwischenbedingungen zu beriicksichtigen, was wieder 
auf Gleichungen der Form (22) bzw. (24) fiihrt, nur daf die in ihnen auftretenden Matrizen ‘3 und Q 
aus mehr als zwei Zeilen bestehen. Ihre Elimination ist indessen nicht wesentlich umstandlicher als 
bei nur zwei Unbekannten, da die Gleichungsmatrix ‘3 
infolge stufenweisen Kintretens der neuen Unbekann- 
ten von ,,fastgestaffelter “Form ist, d. h. oberhalb 
i=0 7 e gi der Hauptdiagonalen nur eine Reihe mit von Null 
verschiedenen Elementen aufweist. Die Berechnung 
der Matrix 2{ verlauft in der gleichen Weise wie oben. 

Wir erlautern das Gesagte an einem Beispiel einer Welle mit Zwischenstiitze und Zwischen- 
gelenk (Abb. 4). 

Die Rechnung erfolgt nach dem Schema 6, wo die — in sich gleichen — Belastungsmatrizen 
fortgelassen sind. 


Abb. 4. Welle mit Zwischenstiitze und Zwischengelenk. 


Rechenschema 6 


Po Q 4Q Age A po A yy Aw, Ao A Age 
ea | 
1 0 
Fi 0 0 Po 
0 1 
Deel soenl 1 1 0 10, ees 0 
(enGrs 1 3 0 Foe ni A 
ran 0 1 0 Ate Sg 
1 0 1 1 Sy SR 
Pele Sch 2 8 1 0 UR) hy ya ib, 
feos 1 12 3 (1) 213 8168 LS all Po, APs 
tect 0 2 1 0 56) =126 6304 0 ; 
1 0 1 1 0 35 0 210 ee =a5 
[ee hes ae) 3 27 8 1 626 —344 585 == .92 10 0 
rou 1 27 12 l 714 =—336 1428 —168 42 0 
Pal 0 3 2 0 9h Ls 28 490 9 oe 21 
(tbs 0 1 1 0 | 35 0 4.20 0 a5 
Endglei-° | 1 1 0 Qe Lom eae 0 0 0 Bee 
chungen 0 2 1 Oi MQ == Ae 63° — 14 0 My 
3 Da 8 1 626 —344 558 oe 10 w3 = 
1 27 12 boss] ee EEL 1428 —168 42 Q3 = 
Elimsnation +) “Wey 0 pe tea Onn tl Coo el estos meaner 0 Pamala Ni sa faa 
Ot 1 0 56) 28 oa haa 0 
3 12, — 4 1 EO a 168 16 10 
i 13 1/4) 3/4 | — 5 31 6s) 42> 5 Ses 
| 28/3 —58/3 — 52/3 20/3 1 
| 44/3 56/3 — 28/3 —124/3 1 
c | —98 98 —259 —868 1 
| 10/3 —10/3 - 62/3 20/3 l 
—16/3 16/3 — 32/3 —158/3 1 
28/3 —44/3 —98 10/3 —16/3 | 
20/355 161/38 196 —20/3 32/3 
% | —10/3 ~56/3 344 —10/3 64/3 


| — 5/3 40/8 385 13/3 80/3 
20/3 —124/3 —868 20/3 —158/3 | 


_ Hier ist als fiinfte BiegegréBe in 1 die Winkelanderung Aq, Ae Stelle des Sees Winkels 
Po + Ap, gewahlt worden, 
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Die Iteration an der Matrix 3 2 zeigt das Rechenschema 7. 


Rechenschema 7 


407 


71 Ws P2 Ags 
ee 
28 — 44 — 294 10 16 
— 20 61 588 — 20 32 
— 10 56 1032 — 10 64 
— $5 40 1155 13 80 
20 —]24 —2604 20 — 158 
— 1 1 1 1 2 
—324 585 960 1053 —2412 
— 0,3 | 0,6 1 i 2.5 
—278,8 | 530.6 898,60 993,5 —2209,4 
— 0,31 0,59 il 1,11 200 
—277,06 927,03 895,12 992,18 —2261,42 
— 0,309 523 0,588 781 | 1 1,108 432 2,526 388 
—277,066 | 927,093 | 895,294 992,397 —2261,861 
— 0,309 469 | 0,588 737 | i 1,108 459 2,026 389 
—277,063 527,089 | 895,290 992,396 —2261,854 
— 0,309 467 0,588 735 il 1,108 463 2,020,393 — 1D, 
895,290 
32, = 895,290 
= 298,430 
At = 8,444 19 
A, = 1,704 67 
exakt 1,701 50 
Fehler 0,003 17 = 1,9°/o9 


In Anbetracht der groben Unterteilung und des bescheidenen Rechenaufwandes diirfen die 
Ergebnisse als sehr gut bezeichnet werden. 


(Eingegangen am 17, Januar 1958.) 


Anschrift des Verfassers: Dr.-Ing. R. Zurmuhl, Darmstadt, Martinstr, 75. 
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Some Special Solutions for the Flexural Vibration 


of Dises of Varying Thickness 
By H. D. Conway 


1. Object of Investigation. The present investigation concerns the calculation of the frequencies 
of axisymmetrical vibration of discs whose flexural rigidity D varies with the radius r according to 
the law D = D)r", where Dy and m are constants. Although it was not found possible to easily 
solve the problem for arbitrary m, some particular solutions in terms of Bessel functions were found 
for a number of positive values of m and corresponding values of Poisson’s ratio y. Although these 
solutions have limited practical value, they are nevertheless exact and hence can be used to assess 
the accuracy of results obtainable by approximate methods. Numerical results are obtained for 
plates having a clamped edge. Other investigations relating to this topic have previously been 
given by the author” ®. 


2. Development of differential equation. The differential equation to be solved is*® 


a /®w 1 dw\ , @D/ew  v» ew 1 f . ew 
: } =a | } es =< Ss 
p or & iS ore Or ' dr (a8 Ue iF | on az” dr, (1) 


where w is the deflection at radius r, D is the flexural rigidity, h is the thickness, 9 is the density and 
y is Poisson’s ratio. 

In order to obtain equations which are tractable, we tentatively write the flexural rigidity D 
as proportional to a power of the radius r and hence D = Dy r™, where Dy = E hi/12 (1 — »”) and 
m are constants. Thus the thickness h of the disc is given by h = h, r/3. Equation (1) then becomes 


o*w ; ay CEH 0?w Ow 
ie a r? (2 m + 2) ae +r? (m + ym —1 + m?) 78 +r (1 + »m*—vm—m) ir 
12 9 (1 — 2”) (12 —2 m)/3 Pw : 
cob hte S de (2) 


0 
For vibrations, we further write 


w = f(r) er (3) 


and hence obtain the following equation for f(r): 


4 d3 r d df 
Aci trem +2) +2 (m4 ym— 1 + m?) Lt +m? —ym—m) 7 


op (Lr) ay oes 
ee r(2—2m)/3 f — 0), (4) 


The solution of this equation will now be examined. 


3. Solution of differential equation. The general solution of Equation (4) is obviously difficult 
to obtain. However, consider the following equation given in the book by Kamke* 
dif 
4 


df a d 
mt Ars t+ 44+ 47L£4 Af=0, (5) 


r 


where 

A, =6—4a—4c, 

A, = 2 (a? — pc?) + 4 (a + ¢— 1)? + 4(a— 1) (e—1)— 1, 
A, = [2 (wc? — a?) — (2a— 1) (2e—)] (2 a+2c—1), 
Ay = (a — wc?) (a? + 4a + 4c? — 2 c2) — Bt ct re, 


Its solution is 
a =r [C, J(u) + Cy Y,(u) + C, J,Aiu) + Gy Y,(-u)], wath we bro (6) 

* H. D. Conway, J. Appl. Mech., June (1955) p. 275. 

; a D. Conway, Z. angew. Math. Mech, 37 (1957) p. 406. 

“E, 


Timoshenko, Theory of Plates and Shells, p. 283, New York 1940, 
Kamke, Differentialgleichungen, p, 534, Leipzig 1943. 
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Writing a = we = — m/3, c = 1 — m/6, wp = 2 m|/(m — 6), Equation (5) becomes 
ee + 73 (2m + 2) = + r( ae L ah ) is 
as r(t wit ae m) ws b(1 +) p12 —2m)/3 f — 0), (7) 
Further, writing Poisson’s ratio y = (2 m th 3)/9 in RaWition (4), we obtain 
iT +8 (2m 42) oF 4 32 Ges ae a 
risa HEP moyn, 


Equations (7) and (8) are identical if b4 (1 — m/6)! = 120 p? mh —v)/E h? and v = (2 m— 3)/9. 
Considerable restriction is thus placed on the solution when it is eomcushered that » must lie in the 
range 0 — 1/2. However, the results will show the effect of certain tapers on the frequencies of 
vibration and, since they are exact, they may be used to test the accuracy of results obtained from 
approximate methods developed to investigate the vibration of more practical discs. The following 
special cases are now considered. 


4. Special cases. a) m = 2, y = 1/9. The solution of Equation (4) is here 
OYE , - ‘ he 5 i 600 p* 
J=r *(C, J_u) +¢,Y_ vy) + C, Juitu) + ¢C,Y_x(iu)], with uw=br" and t= pee 
0 
For a solid disc we note that the second and fourth terms in parentheses are singular at r = 0. 
Hence C, = C, = 0 and, since J_,(u) = J,(u), J_i(i u) = i I,(u), we write with no loss of generality 


f=r*8[B, Jy(u) + By T(u)]. (10) 
For a solid disc clamped at its outer radius r = ro, it follows that the equation for the determination 
of the circular frequencies p is 
Jiu) Liu) — bri, 1s okey ll 
tian ha with uy = bré and b Eke’ (11) 
the primes denoting differentiation with respect to r. It is interesting to note that this is also the 
equation the finding of the first root of which is necessary to determine the frequency of a uniform 
thickness disc vibrating with one nodal diameter and no nodal circles, that is, the second harmonic 
of the disc. 
The first ten roots of Equation (11) have been tabulated by Airey}, the first root being uj = 450.4. 
The fundamental pericd for the varying-thickness disc is then found and, denoting the thickness at 
the outer edge by H, is given by 


- (9) 


60 9 p* 


dx _ Baa 1 yo (12 
bi, eH 8579 | 

The corresponding result for a constant thickness disc has, for y = 1/9, factor of 38.74 instead 

of 35.79, and it follows that the ratio of the fundamental period of the variable (D = Dr’) 

to the constant thickness, clamped-edge disc is 1.08. 


b) m = 18/7, vy = 5/21. The solution for this case is 
= 637 
je roe [G J_3/2 (u) aL C, Y_3/2 (u) = OF J_3)2 (i u) 3 C,Y_3)2 (i u)| , u=b ie Us 6 a (13) 


For a solid disc, the first and third terms in parentheses are singular at the center, and noting that 


Y_3/2 (u) = ee eee + cos u). Yow = 1 oe +- cosh u (14) 


tu u TU uw 


the expression for f may be written in closed form. For a clamped disc of radius rp this leads to 
the frequency equation 


> We, i 637 @ p? y 
u, (tan uy + tanhu,) =2tanu,tanhu,, with u,=—Obril? and b= 6 ER (15) 


The first root of this equation is uy = 5.27, the fundamental period is found and, denoting the 
thickness at the outer radius by H, is given by 


2% ~— 32207 1 ; a, (16) 
Pp He ss21-7 3k 


44, JLR Re Airey, Proc. Phys. Soc. London 23 (1911) pp. 227, 231, 232. 
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The corresponding result for a constant thickness disc has, for »y = 5/21, a factor of 39.65 instead 
of 35.21, and thus the ratio of the fundamental period of the variable (_D = D, r'®/") to the 
constant thickness, clamped-edge disc is 1.13. 

c) m=3,y =1/3. This particular case has been discussed before by the author! as presenting 
an interesting analogy between the symmetrical vibrations of a disc of linearly varying thickness 
and a cone. We quote the results here. The solution is 
far [C, J_a(u) + GY_2(u) + C,J_o(iu) +C,¥_2(iw)], with u=br? and b= shH (17) 


The equation for the determination of the frequencies is 


J(u) Ly(uo) = Jato) To(uo) with uy = bri’. (18) 
The first root is uy = 5.906 and the fundamental period is 
2 32%a* 1 jee (19) 
P H 34.88 V3£E 


The corresponding result for a constant thickness disc has, for y = 1/3, a factor of 40.84 instead 
of 34.88, and thus the ratio of the fundamental period of the variable (D = Dy r*) to the constant 
thickness, clamped-edge disc is 1.17. 

d) m = 6, y arbitrary. The writing of m = 6 in Equation (4) makes the latter homogeneous, 
and with the transformation r = e*, the equation becomes 


Se 2 Rea) See ee reesiye (20) 
It is then necessary to find the roots of the quartic characteristic equation 
4828 +2(5 +39) 2—24(1—v)2Z-—K=0 with K= aa a (21) 
The above equation may be written in the form 
(A +2 + B,) (A + 2—f,) (A + 2 +B) (A + 2 — fy) = 0 (22) 
whence, by multiplication and comparison with Equation (21), it may be shown that 
+A =M—6v, (fff? = 36(1—v)? 44 K. (23) 


Thus the roots of the characteristic equation are 


v42=—2ty7—3v+/9— 9 EK, yn4=—24V7—39-f9— PF EK (24) 


and the equation for f is 


f=Arm+ Bre+Cre+ Dr. (25) 

However, it will be noted that since z = log r, which is singular at r = 0, the above solution is 

valid only for a disc with a hole. If the vibrations of a disc with a clamped outer and free inner 

edge are investigated, it is found, as might be imagined, that the frequency decreases indefinitely 
with decrease in the radius of the hole. 

5. Summary. The free, flexural, axisymmetrical vibrations of a clamped circular disc having 

a flexural rigidity given by D = Dy r” have been investigated, solutions in terms of Bessel functions 

being found for certain values of m and Poisson’s ratio y. Values of the ratio of the fundamental 

period T of the varying thickness discs to the fundamental period T, of corresponding constant 


thickness discs having the same value of Poisson’s ratio have been computed, and are tabulated 
below 


| 
m | v GD 
2 Gow ic0s 
245 |. “5/2t | 1.48 
Shh ls acer ig 
The case of m = 6 (parabolic variation of thickness) and arbitrary v has also been investigated. 


(Eingegangen am 25. Januar 1958.) 


Anschrift des Verfassers: Prof. H. D. Conway, Ithaca, N. Y. (USA), Cornell University. 


' Vide reference in footnote 2 on page 408. 
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Uberschallstrémung am Austritt von Turbinenschaufel-Krinzen 


Von H. Hausenblas 


1. Einleitung. Wie Stodolas Ubersicht! zeigt, wurden die Dampfturbinenbauer schon bald zum 
Studium der Beschleunigung einer Strémung auf Uberschallgeschwindigkeiten im Austrittsdreieck 
von Turbinenbeschaufelungen angeregt. Eine kurze, ziemlich vollzahlige Zusammenstellung der 
verschiedenen Arbeiten gab kiirzlich Linnecken®. Hierin nicht enthalten ist eine neuere Arbeit 


von Traupel*, die einen wesentlichen Beitrag zum Verstandnis der eindimensionalen Rechnung 
dieses Vorganges enthalt. 


SS 
oe ; SG, Austrittsenden der Schauteln 


yy \ Ib Verdlinnungslinien 
Tee 


poe —+ ct ie o a ae ie te Stromlinien 
) Ss | : eo ae) 
Ne x —% \ Y. SS & l \ " 


ee oS N 
fee Che he oe ee 
Verdichtungsstok Wirbelfldchen 


Abb. 1. Uberschallstrémung hinter einem Streckenprofilgitter. 


Hier soll das Problem nur fiir den Fall erértert werden, daB® die (Uberschall-)Strémung auch 
hinter dem Austrittsdreieck der Beschaufelung durch zwei koaxiale Zylinder begrenzt wird (innen 
Nabe, auBen Turbinengehause*). Nach der tiblichen Abwicklung des betrachteten Zylinderschnittes 
der vollbeaufschlagten Axialbeschaufelung in ein gerades Schaufelgitter ist also ebene Strémung 
vorauszusetzen. ‘ 

Bei der Auslegung von Turbinen interessiert man sich praktisch nur fiir die auftretenden 
Austrittsablenkungen. Einfache und bequeme Berechnungsmethoden hierfiir erhalt man aus einer 
eindimensionalen Behandlung der Uberschallstrémung im dreiecksférmigen Austrittsabschnitt 
eines Beschleunigungsgitters. Fiir das nahere Verstandnis dieser Strémungsvorginge ist aber eine 


A, Stodola, Dampf- und Gasturbinen, Berlin 1924. 
H. Linnecken, AEG-Mitt. 47 (1957) S. 44. 


1 

3 W. Traupel, Die Strahlablenkung in der vollbeaufschlagten Turbine, Mitt. Inst. f. Therm, Turbomasch. 
ETH Ziirich Nr. 3, Ziirich 1956. 
4 


H. Hausenblas, Energie 10 (1958) S. 131. 
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zweidimensionale Betrachtung erforderlich. Die Darstellung von Fliigel! ist nicht ganz korrekt, 
wie ich zeigen konnte. Auch Meyer® gibt mehr anschauliche als exakt begriindete Darlegungen. 


2, Anwendung des Charakteristikenverfahrens. Eine exakte zweidimensionale Behandlung des 
Uberschall-Strémungsfeldes im Austrittsdreieck ist fiir ebene, reibungsfreie Stromungen mittels des 
Prandtl- Busemannschen Charakteristikenverfahrens méglich®. Abb. 1 zeigt nochmals das von mir 
erhaltene Ergebnis. Es wurde beriicksichtigt, daB sich die aus den einzelnen Schaufelkandlen 
kommenden Strémungen hinter den Austrittsdreiecken zu einem geschlossenen ebenen Stro- 
mungsfeld vereinigen, diese sich also gegenseitig beeinflussen. Damit die Strémung dem ebenen 
Schaufelriicken folgen kann, werden die von der Schaufelhinterkante A ausgehenden Verdiinnungs- 
linien am Schaufelriicken der Nachbarschaufel reflektiert und laufen aus dem Austrittsdreieck 
stromabwirts heraus. Diese Reflexion bedingt den Unterschied der Strémung im Austrittsdreieck 
von der Strémung um eine konvexe Ecke nach Prandtl-Meyer*. Langs des Schaufelriickens tritt 
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Abb, 2. Uberschallstrémung hinter einem speziellen Gitteraus gekriimmten dicken Profilen mit scharfer Austrittskante. 


eine starkere Expansion ein. Da hierdurch die Strémungen zu beiden Seiten der Schaufelhinter- 
kante verschiedene Machzahlen und verschiedene Abstrémrichtungen aufweisen, mi ten sie 
einander durchdringen. Dies wird vermieden, indem von der Sehautcihineerkenes ein Verdich- 
tungsstoB ausgeht, der die Uberexpansion der Strémung langs des Schaufelriickens (verlustbe- 
haftet) riickgangig macht®. Dieser VerdichtungsstoB hért aber nicht, wie Meyer vermutet, mitten 
im Strémungsfeld auf, sondern lauft mit (durch die Einmiindung der Verdannnnethnen ab- 
nehmender Starke stromab. Die Abstrémung des Gitters durchquert die von den verschiedenen 
Schaufelhinterkanten ausgehenden VerdichtungsstéBe, wobei die StoBstarke von Kreuzungsstelle 


1G. Fliigel, Die Dampfturbinen, S. 67, Leipzig 1931 
2 H. Meyer, Z. VDI97 (1955) S294, 
: ae FuBnote 4 von S. 411. 
. Sauer, Theoretische Einfuhrung in die Gasd i 1 
* R. . ynamik, S. 58 und S. 108, B 
* Die Isentropie der Stromung endet natiirlich an VerdichtungsstéBen. ae 
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zu Kreuzungsstelle abnimmt. Dieses Ergebnis ist insoweit noch unbefriedigend, als ein Strecken- 
profilgitter vorausgesetzt werden muBte, dessen ,,Schaufelkanale“ mit Schallgeschwindigkeit durch- 
strémt werden. Bei der Wahl anderer Schaufelformen mit beschleunigter Strémung liegt die 
Schwierigkeit in der Berechnung der schallnahen Unterschallstrémung und der Festlegung der 
Linie, auf der die Schallgeschwindigkeit durchschritten wird (Schall-Linie). 

Um wenigsiens ein Beispiel fiir eine allgemeinere Schaufelkanalform zu erhalten, muB man 
bis zur Schall-Linie einen Ausschnitt aus einer bekannten beschleunigten Strémung mit Durch- 
gang durch die Schallgeschwindigkeit beniitzen. Fiir Abb. 2 wurde eine Wirbelsenkenstrémung 
(Mittelpunkt 0) gewahlt, die durch Linearverbindung von Senkenstrémung und Wirbelstrémung 


Senke Wirbelsenke Wirbe/ 


Wirbelsenke 


Abb. 3. Erzeugung einer Wirbelsenke durch Linearyerbindung einer Senken- und einer Wirbelstrémung (Teilverhaltnis = 1). 


entsteht. Einfachheitshalber wurde das Teilverhaltnis! ,/A, gleich Eins genommen. Sind r, 
und r,, die Halbmesser fiir Kreise derselben Machzahl M* der Senken- baw. Wirbelstrémung, so 
erhalt man nach Abb. 3a: 


iy Gaara e 
Boe ne ae tg) =—, 


wobei r der Halbmesser des Kreises der Wirbelsenkenstrémung, auf dem die Machzahl M* herrscht, 
? die Neigung des Geschwindigkeitsvektors gegeniiber der Umfangsrichtung des genannten Kreises, 
M* die Machzahl bezogen auf die kritische Schallgeschwindigkeit, 0 die zu M* gehorige 
Stromdichte, dimensionslos gemacht durch Bezugnahme auf ihren kritischen Wert, ist und 
C, sowie C, zwei Konstanten bedeuten. Fiir den Zentriwinkel p der spiraligen Stromlinien ergibt 


sich (Abb. 3b): 
eld " dr Z 
ee ripe’ 


"a 


der Index a kennzeichnet irgend einen Anfangswert. 
Fiir Abb. 22 wurden die Kreise von Senke und Wirbel fiir M* = 1 gleich groB gewahlt, so dab 
die Stromlinien der Wirbelsenke diesen Kreis unter 45° schneiden. Der Zentriwinkel des beniitzten 


1 Siehe FuBnote 4 von S. 412. 
2 In Abb. 2 wurden gegeniiber Abb, 1 die Wirbelschichten nicht mehr vermerkt. 
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Ausschnittes der Wirbelsenkenstrémung (16°) konnte nicht beliebig groB gewahlt werden. Es 
wurde stillschweigend angenommen, da die Strémungszustande weiter stromauf (z. B. in der 
Umgebung der Schaufeleintrittskanten), die von der Wirbelsenkenstrémung abweichen, letztere 
nicht beeinflussen, was bei engen Teilungen ausreichend zutrifft. Obwohl sich die berechnete 
Wirbelsenkenstrémung bis zu einem Grenzradius in den Uberschallbereich erstreckt, wurden die 
numerischen Ergebnisse nur fiir den entsprechenden Ausschnitt des Unterschallbereiches tiber- 
nommen. Fir die Fortsetzung der Strémung von der Schall-Linie in den Uberschallbereich diente 
das Prandil-Busemannsche Charakteristikenverfahren. Dabei wurde vom Uberschallbereich der 
Wirbelsenkenstrémung nur eine dreiecksférmiger Teil zwischen der kreisférmigen Schall-Linie 
ABCDE und dem Linienzug AFGHJ beniitzt. Die Schall-Linie wurde von der Schaufelhinter- 
kante A ausgehend angenommen, da sonst ein mit Schallgeschwindigkeit durchstrémter parallel- 
wandiger Teil des Schaufelkanals aufgetreten ware. 

Da hier wegen der geschilderten Schwierigkeiten die Schall-Linie von vornherein festgelegt 
wurde, war damit die Form des Schaufelriickens auch stromab vom Punkt EF weitestgehend fest- 
gelegt. Der Bereich von E bis J muB der Stromlinie der Wirbelsenkenstrémung entsprechen. Um 
in den Punkten B, C und D die fiir die Wirbelsenkenstr6mung notwendigen Verdichtungslinien 
in die Schall-Linie einmiinden zu lassen, mu® der Schaufelriicken die Knicke K, ZL und M er- 
halten. Von der Anbringung weiterer Knicke in diesem Bereich wurde abgesehen, obwohl dies 
méglich erscheint. Von diesen wiirden stromab entsprechende Machsche Linien ausgehen, die 
auch die Strémung hinter dem Gitter beeinflussen wiirden. Von M an wird in Ubereinstimmung 
mit Abb. 1 der Schaufelriicken geradlinig angenommen. Die dadurch bedingten Vorgange in der 
Strémung stimmen vollkommen mit den fiir Abb. 1 friiher von mir beschriebenen iiberein?. 

Die in Abb. 1 und 2 gezeigten Strémungen entsprechen jeweils dem niedrigsten Gegendruck, 
der durch Expansion im Austrittsdreieck erreichbar ist. Die letzte von der Schaufelhinterkante 
ausgehende Machsche Linie trifft dabei gerade auf die Hinterkante der Nachbarschaufel. Bei 
hdheren Gegendriicken verbindet die letzte Verdiinnungslinie im Austrittsdreieck nicht mehr die 
Austrittskanten benachbarter Schaufeln, sondern miindet auf dem Schaufelriicken weiter strom- 
auf. Der von der Austrittskante ausgehende VerdichtungsstoB wird entsprechend schwacher und 
verschwindet bei Abstrémung mit Schallgeschwindigkeit. 


3. Allgemeine Ergebnisse. Das untersuchte spezielle Beispiel 1aBt die folgenden allgemeinen 
Ergebnisse erkennen: 

Am gekriimmten Schaufelriicken von Abb. 2 ergeben sich im Gegensatz zum geradlinigen 
von Abb. 1 im Gebiet gegeniiber der benachbarten Austrittskante 4 bekannterweise Uberge- 
schwindigkeiten. Der Schalldurchgang findet daher am Riicken nicht mehr querab zur Austritts- 
kante A statt, sondern weiter stromauf. Die besondere Form der Schall-Linie (hier ein Kreis- 
bogen) und der Strémung im Gebiet 4EM hangt von den besonderen Annahmen unseres Bei- 
spieles ab, sind also nicht allgemein giiltig. Hierzu gehért auch die Tatsache, daB die Krimmung 
des Schaufelriickens in der Umgebung von J am gréBten ist. Die in Abb. 1 gezeigte Strémung 
im Austrittsdreieck findet sich in Abb. 2 unverandert wieder. Dies ware allerdings nicht mehr 
der Fall, wenn, wie oben erwahnt, die Kriimmung des Schaufelriickens im Bereich J bis M anders 
gewahlt worden ware, so dafS von diesem Bereich zusatzliche Machsche Linien ausgehen, die dann 
auch die Strémung hinter dem Gitter beeinflussen wiirden. 

Wie schon erwahnt, wiirde im Normalfall des hier behandelten Problems nicht die Form der 
Schall-Linie, sondern die Gestalt der Schaufelkontur (vor allem des Schaufelriickens) vorgegeben 
sein. Die Bestimmung der zugehérigen Strémung (insbes. der Schall-Linie) sté®t aber auf erheb- 
liche Schwierigkeiten. 

Das Beispiel Abb. 2 enthalt ferner die Beschrankung auf eine scharfe Hinterkante. Bei rei- 
bungsfreiem Gas ware hinter einer dicken Austrittskante (z.B. mit halbkreisférmigem Ende) 
eine gemischte Uberschall-Unterschall-Strémung zu erwarten. Die Strémung wiirde sich bei 
Umstrémung des konvexen Teiles der Austrittskante auf héhere (Uberschall-)Geschwindigkeiten 
beschleunigen, um sich dann in zwei schrag ablaufenden VerdichtungsstéSen wieder zu verlang- 
samen. Der dem Austrittsdreieck abgewandte VerdichtungsstoB wiirde sich wahrscheinlich mit 
dem von der Schaufelaustrittskante ausgehenden VerdichtungsstoB der Abb. 1 und 2 zu einem 
einzigen StoB vereinigen. Kine genaue Untersuchung dieser Strémung bereitet wegen ihres trans- 
sonischen Charakters Schwierigkeiten und lohnt nicht, da sich bei den wirklichen, reibungsbehafteten 
Gasen die Strémung von einer solchen dicken Hinterkante ablist. Es entsteht ein Nachlaufgebiet. 


1 Siehe FuBnote 4 von Seite 411, 
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Dies erkennt man deutlich in den Schlierenaufnahmen! von Luksch und Quick?. Dort sieht man 
auch, daf} trotz der Ausbildung einer Nachlaufzone, wie oben beschrieben von den dicken Schaufel- 
austrittskanten jeweils zwei VerdichtungsstéBe ausgehen. Dies la8t vermuten, daB die Nachlauf- 
zonen kurz hinter den Schaufelaustrittskanten schmialer als diese sind. 

_ Eine fiir die Anwendung in der Praxis ausreichend einfache Berechnung des vergleichmaBigten 
Uberschall-Strémungszustandes weit hinter dem Gitter auf Grund der hier gezeigten Unter- 
suchungen diirfte wohl kaum médglich sein. Dort wird man sich also nach wie vor mit den ein- 
dimensionalen Rechnungen begniigen miissen. Wie die Versuchsergebnisse von Wewerka und 
Schmidt? zeigen, geniigen diese in der Regel den von der Praxis gestellten Anforderungen. Aller- 
dings bestehen Schwierigkeiten in der Berechnung des niedrigsten, durch Expansion im Austritts- 
dreieck erreichbaren Gegendruckes, wie ein Vergleich mit den diesbeziiglichen Darlegungen von 
Traupel* zeigt. Dieser Nachteil wird dadurch abgeschwacht, daS (zumindest in der isentropen 
Naherung) bereits oberhalb dieses Gegendruckes jener Punkt erreicht wird, wo bei Steigerung 
des an die Turbine angelegten isentropen (adiabat.) Gefalles die Wellenleistung wieder abzu- 
nehmen beginnt®. 


4. Zusammenfassung. Es wird die reibungsfreie Uberschallstrémung im Austrittsdreieck eines 
Beschleunigungsgitters fiir den Fall eines Streckenprofilgitters und eines speziellen Gitters aus 
gekriimmten, dicken Profilen, jedoch mit scharfer Austrittskante, mit Hilfe des Charakteristiken- 
verfahrens ermittelt. Der Fall einer dicken Austrittskante und die eindimensionalen Naherungs- 
rechnungen werden kurz diskutiert. 


(Eingegangen am 22. Februar 1958.) 


Anschrift des Verfassers: Privatdozent Dr.-Ing. H. Hausenblas, Dudweiler (Saar) Gartnerstr. 2a. 


1 Diese Aufnahmen wurden auf dem vom Verfasser erstellten Uberschall-Gitterwindkanal der SNECMA/ 
Paris erhalten. Bei der Auswertung dieser Schlierenaufnahmen ist zu beachten, dafs meist unterhalb des frei- 
Jiegenden Schaufelkanals nur mehr ein weiterer folgt. An der Strahlgrenze dieses Schaufelkanals werden bei 
diesen Versuchen die Verdichtungsstofe reflektiert, was im ringformiegn zusammenhangenden Stromungs- 
kanal einer Turbomaschine nicht auftritt. ; 

2 W. Luksch u. A. W. Quick, Bullet. Assoc. Techn. Marit. et Aéron., Session 1953, 

3 Wewerka u, Schmidt, Untersuchung von Lavaldiisen und Mindungen; Forschungsber. 1409 d. Deutsch. 
Luftfahrtforschg. 

4 Siehe FuBnote 3 von Seite 411. 

5 Siehe FuBnote 4 von Seite 411. 
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Stationire Warmespannungen mit temperaturabhangigen Stoffwerten 
Von R. Trostel 


1. Allgemeines. Als Ergainzung meiner fritheren Arbeit+, in der das Problem der Berechnung 
stationarer Warmespannungen im mittleren Bereich eines langen Hohlzylinders aus elastisch- 
inkompressiblem Material untersucht worden ist, wird nachfolgend der Weg einer allgemeinen 
Lésung stationirer Warmespannungsprobleme bei temperaturabhangigen Stoffwerten FE, x, und A 
und fir beliebige, jedoch mit der Temperatur nicht veranderliche Querkontraktionswerte y auf- 
gezeigt. Ausgehend von der allgemeinen thermisch-elastischen Verschiebungsgleichung wird, 
nachdem man den Verschiebungsvektor nach einem, die Temperaturveradnderlichkeit des Elastizi- 
tatsmoduls E enthaltenden Parameter entwickelt hat, ein System von Differentialgleichungen fiir 
die einzelnen Glieder der den Verschiebungsvektor darstellenden Funktionenreihe gefolgert. Dieses 
Differentialgleichungssystem ist rekursiv lésbar, sofern das Temperaturfeld und die Temperatur- 
abhingigkeiten EF = E(9) und a, =a,(0) des Elastizitaétsmoduls und des Warmeausdehnungs- 
koeffizienten vorgegeben sind. Zur Ermittlung des Temperaturfeldes gehen wir vom Fourierschen 
Theorem aus, und erhalten aus der Warmemengenbilanz am Kérperelement bei Beachtung der 
Temperaturveranderlichkeit der Warmeleitzahl/ eine nichtlineare Differentialgleichung fiir das 
Temperaturfeld 9, die jedoch (fiir beliebige A(9)) vermége einer geeigneten Transformation auf 
eine lineare Differentialgleichung zuriickgefiihrt werden kann. 

Im AnschluB an diese in Ziff. 2 und 3 dargestellten Untersuchungen folgen in Ziff. 4 zwei Bei- 
spiele, wovon vor allem das zweite, die Stabbiegung bei gleichzeitiger Erwarmung, von Bedeutung 
ist. Man erkennt hier, daB man auch im Falle der Temperaturabhangigkeit des Elastizitatsmoduls 
die Bernoullische Hypothese vom Ebenbleiben der Querschnitte aufrechterhalten kann, und dah 
in jedem statisch bestimmt gelagerten Stabe bei ungleichmaRiger Krwarmung Temperaturspan- 
nungen entstehen. Sie sind jedoch in der Regel geringfiigig. 

Das andere Beispiel behandelt wieder den Hohlzylinder. Hier bestatigen sich die Ergebnisse, 
die der Verfasser in dem bereits erwahnten Aufsatz vorausgesagt hatte: Die tatsachlichen Warme- 
spannungen in einem Stahlrohr sind gréfer als diejenigen, die sich unter der Voraussetzung un- 
veranderlicher, bei Zimmertemperatur vorhandener Materialwerte ergeben. Fiir die tatsachliche 
Warmespannungsverteilung wurden hinreichend genaue Naherungsformeln hergeleitet. Mit ihrer 
Hilfe erkennt man auBerdem, daB man den EinfluB der Temperaturveranderlichkeit der Stoff- 
werte EF und «, ziemlich genau auch mittels der Theorie unveranderlicher, fiir eine mittlere 
Systemtemperatur #,, giiltiger Materialwerte E,,, und «,,, erfassen kann, sofern die Temperatur- 
unterschiede innerhalb des Systems nicht extrem grof sind. 


2. Die Ermittlung des Spannungs- und Verschiebungszustandes. Wir setzen zunachst eine be- 
hebige Temperaturveranderlichkeit des Elastizitatsmoduls EF und der Warmeausdehnungszahl «, 
voraus. Die Querdehnungy wird nicht als temperaturabhangig angesehen. Die maBgebenden 
Gleichungen sind die Gleichgewichtsbedingungen fiir das Kérperelement, die mit dem Spannungs- 
tensor © und den auf die Volumeneinheit bezogenen Massenkraften f die Form 

VY OA =0 (1) 
annehmen, sowie die durch die Temperaturdehnungen erweiterten Hookeschen Gesetze, die man 


in der Form 
0 


x v : 1l+y : J 
See se { a,(t) dt) & (2) 


6 


schreiben kann, wobei mit dem Verschiebungsvektor 
Lie a - 
Rete ees) (3 


der Verzerrungstensor?, © der Einheitstensor und # die Temperatur gegeniiber dem spannungs- 
losen Ausgangszustand sind. 


 R. Trostel, Ing.-Arch. 26 (1958) S, 134. 
* Das Zeichen bedeutet die dyadische Multiplikation. 
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Einsetzen von (2) in (1) liefert die Verschiebungsgleichung des Problems: 


[ 0 
y ; lty . 
V e (57 (ee Oy a6 


0 


+4, (VE)|D 


a 
| [ees 7 : 
(7 — | OX, a) i (4) 


Nun ist 


VD=>FVIVou+ucoy] =i CY) yl = Aut Vv), 


ViCy ay el = 7 ws = V¥ (Wu), V[(fow)e]=v soa, 


und setzt man 


ee mF E) = « Bid), (5a) 


E do dé 


wobei € ein die Konstanten des Materialgesetzes F = E() enthaltender Parameter ist, so daB 
sich also 


1 1 dE 
VE =e ZVI =e) VI=cr, r= OVE (5b) 
ergibt, so folgt schlieBlich aus (4) die Relation! 
+ \ v 
oa) Hei + Gy +E ea 7 
a 
2(1 +» 2 é 
eo (6) 


die fiir ¢ = 0 sowie o, = 0,9 in die thermisch-elastische Verschiebungsgleichung nach der Theorie 
unveranderlicher Materialwerte tibergeht. Denkt man sich die Lésung von (6) im Sinne der Me- 
thode von Poincare in Form der Entwicklung 


u= Yu, =U t+ > ey, (7) 
k=0 hat 


dargestellt, so liefert Einsetzen von (7) in (6) nach Koeffizientenvergleich hinsichtlich e* das 
Gleichungssystem 


Aut rags V (td = 30599 fo ae — 2, (8a) 


Ai -q— 55 V (Vth) =— lambs elor, SG = 10 fa sal (8b) 
L 0 


At, +=, V (VW) =—| VJ) ite + (Ete) 52 (0 Ta) (k=2...n). (8c) 


Mit (7) nehmen nun nach (2) die gee 2 die Gestalt 


ae a ae rand] 
z 0 
ee ee ee b (9) 
k=1 es 
E pee 0 y Fae ; x 
eel (57 ce =i) | ta (VIE). 
4 . 


1 Die Pfeile deuten an, da die Nabla-Operation allein auf den Verschiebungsvektor u anzuwenden ist. 
US 
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an. so daB man hiermit die Gleichungen (8) noch etwas einfacher schreiben kann: 


Ripe a) AY fa ee (10a) 
D 
Ait ea (Vi) =—2(1 +9) 2 = 20 + F(V 9) GS, (10b) 
ne D x 
Ai, +75, V (Vi) = 2.0.4 y) At _ 21 +) S (VO Er-1 (b=2..-m). (100) 


Bei vorgegebenem Temperaturfeld ) = (x, y, 2) sowie bei vorgegebenen Materialgesetzen L = E(?) 
sowie x, =a,(9) ist das Gleichungssystem (8) baw. (10) rekursiv lésbar. Das Randwertproblem 
wird einmalig fiir den Grundzustand (Up, Go) befriedigt. Samtliche iibrigen Randwertprobleme 
fiir die Zusatz-Zustande (U,,©,) baw. (uy, ©,) sind dann homogen. 

In hinreichend diinnen scheibenartigen Gebilden kann ein ebener Spannungszustand 
vorliegen, sofern die Massenkrafte keine Komponente senkrecht zur Scheibenebene hesitzen. 
Sind durch e, und e, zwei in die Scheibenebene fallende orthogonale Richtungen charakterisiert, 
so lat sich der Spannungstensor mit Hilfe der Normal- und Schubspannungen o,, ¢, und T,, des 
ebenen Zustandes in der Form 


ss Onan t 
OS eee a eG eels faa (11a) 
t=, 8 par ip TB x Op 
0 0 
darstellen. Wird als Verzerrungstensor D mit \7 = ¢, Bes rae ee - der Ausdruck 
1 
a= (V site Wen) (11b) 


definiert, wobei 1 = u, ey + Ug eg der in der Scheibenebene liegende Anteil des Verschiebungs- 
vektors ist, so 1aBt sich das Hookesche Gesetz fiir den ebenen Fall anstelle von (2) durch 


o 
eee ie i ae “a c (12) 
I 0 


ausdriicken. Einsetzen von (12) in die Gleichgewichtsbedingung (1) liefert schlieBlich anstelle 
von (6) folgende Verschiebungsgleichung fiir den ebenen Spannungszustand: 


Ce aeee ge Vt ry) — Sb dt 


Oe 
‘ SS ie ce ee geo == (13) 


und mit dem Ansatz (7) folgert man schlieBlich hier nach Koeffizientenvergleich hinsichtlich 
é* anstelle von (8) das Gleichungssystem 


, 1+» 2(1 +7) _ a 2d +), 
Alp TV Wl V [deel (14a) 
j i 
= We) = 1 J 2 ¥ a) 7 
se Ne -e V) Mo + (Ete) V + 58 (V ty) —“ eg { % a. (14b) 
; 0 


{ | 1 5 
\ Uy one ~ V7 (V u,) = --|@ V) Ve—a + (Ete—1) V + ee (V in) (k= 2" ion) (le 


Mit dem Grundspannungszustand 


E a 
So = eer Do ae (r ; VY Uy — — OX, i) & (15 a) 
sowie den Zusatz-Zustanden : 
z= _ E =~ y eRe SE 
or >. oP pee WALD) G| (15b) 
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folgt dann fiir die endgiiltigen os 


S=G oe eG = Gi +5 Sy : (15c) 
und es lassen sich hiermit analog (10) die ates (14) noch in folgender Weise vereinfachen: 
Doig + ee ae VW (V Uy) = —— V ie dis ae ans (16a) 

0 
At, +> V (Vi) =—2 +) 46, =—201 + N2(V HG, (16b) 
A, 4 ana ad ov (Vi) = — 2 (1 + 9) = Ge ae (t +29 Dd) Gp-1 (k=2...n). (16c) 


3. Das Temperaturfeld. Mit dem Fourierschen Theorem 


wonach der je Zeit- und Flacheneinheit vorhandene Warmeflu8 gq [calem~* sec] mit der i. allg. 
temperaturabhangigen Warmeleitzahl A [calcm=! sec! °C] proportional dem Temperaturgra- 
dienten gesetzt werden kann, erhalt man aus der Warmemengenbilanz am Kérperelement im 
stationdren Falle, wenn auBerdem der Volumeneinheit je Zeiteinheit eine Warmemenge Q 
[cal em™? sec] zugefiihrt wird, die Relation 


yYq—0=0', Bhs lite VAY =—Q. (18) 


Fiihrt man hierin anstelle der Temperatur # die Funktion 


iy 


O= A(t) dt mit Ag eA Zeb. dur (t= 00G (19a) 
| 
ein, so dah 
1 
VO=7AIVI (19 b) 
“0 
ist, so folgt aus (18) eine lineare Differentialgleichung zweiter Ordnung fiir 0: 
7 (V0) =AO=—F. (20) 
ho 


Nach Lésung von (20) ist nunmehr bei vorgegebenem Materialgesetz 2 = A() das Temperatur- 
feld } sofort bestimmbar. Setzt man z. B. ein lineares Gesetz in der Form 
MP) = Ay —1, 0 (21a) 


voraus, so hat man nach (19a) 
i 
cee eh ho pee ahi a3 
a | (Ay — Ay 1) dt = 9 — 57, 
t=0 


also eine quadratische Gleichung fiir J, aus der man 


G _hy [Ay \ a, dp 7 Xo Ay a 211 
ame yn) “ia aay hel 2729) om 


findet. Fir 2 =/, = konst. ist nach (19a) O = ¥, so dafs dann (20) die gewohnliche Warme- 
leitungsgleichung fiir konstante Warmeleitzahl darstellt'. 


4. Beispiele. a) Erstes Beispiel. Wir betrachten den mittleren Bereich eines sehr langen 
Stahlrohres mit freien Enden, dessen Mantelflachen r =r; = 12cm baw. r =r, = 18cm die 
pomecreres @; = 350° C baw. 3, = 280° C gegeniiber der Temperatur § = 0° C des spannungs- 


1 Auch aus (21b) laBt sich fiir A(9) = Ay, d.h. A, > 0 die Beziehung O = @ folgern. Fiir A,/A) <1 ist 


Ay ibs Ay ; 4 _ 4o 4 Ay 7a lesa 
yi aa wee also oo, 1—14- ay 
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losen Ausgangszustandes besitzen sollen. Die Temperaturabhangigkeit der Stoffwerte kann im 
Bereich bis etwa 400° © nach Hiitte Bd. I, mit geniigender Genauigkeit durch 


MP) =Ay—A, 9; (22a) 
a(P) = %o + %1?> b20, (22b) 
E(9) = E, — £, 0 — E, ® (22¢) 


approximiert werden, wobei fiir FluBstahl 
Ay =A (8 = 0° C) = 0,12 [cal cm sec °C] , A, = 0,7 - 10 [cal em see °C] , 
Oe oe = 0° C) = Lo Oa a i Gaal Oy 1.0 Oe Cele 
Bp EO 0% G) 12512 108 |Keremes |: B= 2,75 30 | eorensa aaa ss 
EAL [Keene 7G] 

ist. Die (mit guter Naherung konstante) Querkontraktion wird zu y = 1/3 angenommen. 


Wir berechnen zunaichst das Temperaturfeld: Im mittleren Rohrbereich bildet sich ein ebenes 
Temperaturfeld aus, so daB (20) mit Q = 0 sowie 
bee MEAT ee 1rd alan ed 
A= OF r es ( # 


/ 


r dr rdr 
die Form 


_1 d/d0\_ 
AO=zE(F)=° 


annimmt. Durch fortlaufende Integration erhalt man 
@ = G InreC,s 
so daf hiermit aus (21b) endgiiltig fiir das Temperaturfeld 


ho praati ters a 
o=Pl1 \p 2F(Guinr + 69) 


hervorgeht. Die Integrationskonstanten werden aus den Randbedingungen ? = #, fiir r =r; 
und # = #%, fiir r =r, berechnet, und es ergibt sich letztlich 


In lin 
A ae 
Mr) = | = 1 Ay Vi au Ts (1 Ay Oa\? r; 
fe ho } nee ho 2 
ee i 
(23a) 


Eine zahlenmaBige Berechnung zeigt, daB sich das durch 
(23a) gegebene Temperaturfeld nahezu iiberhaupt nicht 
von dem Feld mit konstanter Warmeleitzahl 


ip 


9,n*+ Gln 
8(r) =- Satie i (23 b) 


aes 
aay 


t 
unterscheidet!, so da fiir die nachfolgende Spannungs- 
berechnung das Temperaturfeld nach (23b) zugrundegelegt 
wird. 

Da das Temperaturfeld hier allein von r abhangig ist, 
k6nnen wir einen achsensymmetrischen und ebenfalls allein 
von r abhangigen Spannungs- und Verzerrungszustand 
erwarten, Dariberhinaus bleiben im mittleren Rohr- 
bereich simtliche Rohrquerschnitte, z = konst., wahrend 
der Verformung eben, so da® man fiir die Axialverschie- 
ere bungen w mit den Konstanten 4, und A, 


PP es w = w(z) = 4,2 + A, (24a) 


i“ naa ist i. allg. nur zutreffend, sofern beim Temperaturfeld die Randtemperaturen vorgegeben sind. 
‘8 konnen im Hinblick auf die Warmespannungen nicht zu vernachlassigende Abweichungen dann entstehen, 


sofern fur das it emperaturfeld eine andere Randwertaufgabe vorliegt, wie dies vom Verfasser bereits in seiner 
friiheren Arbeit gezeigt worden ist, 
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folgern kann. Der gesamte Verschiebungszustand eines Punktes im mittleren Zylinderbereich 
ist damit durch den Verschiebungsvektor 


u=ve,+we, =v(r)e, +(4,2-+4,)e, (24b) 


charakterisiert, wenn mit ¢, baw. e, die Einheitsvektoren in Radial- baw. Axialrichtung und mit 
v = v(r) die Radialverschiebungen bezeichnet werden (Abb. 1). Mit 


SO a eee ee 
Ce tor 8 8S, Gp | te 
= 1 oe, 1 
erhalt man dann unter Beachtung von — ag ee 
fs 1 fe 
dv v i daleland 
Yaa e+ +4 =—f +4, Vivw=es E reo, VX =0 


Ata VO V xX (VX) SV a). 
Weiterhin ist 


GQ, ak E, 
la E, +258 2E, 2B, Sse 
Eds”  E,—E,0—E£,& Ey E,o  E, ? pe oe apa 
Ky Ey 
und wir setzen demnach 
o + By ry) ais dp 
ok, x 2 E, a ee 5 Perel Sas) 25 
fae (ee ear erat) pre) 
Ey Ey Ey Eo 
so daB schlieBlich aus (10) mit f = 0 die Gleichungen 
b 
dy(laid 14+ d Ley dé 
ar fra O)| Toya [od =p) S sey 
eid, 1—2»)(1+») Y)a,, (2) eee) 
he a aaa ra Tae eee cay 
. 1—2r)(1+%) YO), ,_ 1—2»)(1 7 
[2d (i)| =< — SSPE) Anat = CEH G3) (k= 2-8) 260) 
hervorgehen. 
Wir betrachten zunachst die Grundlésung. Aus Gleichung (26a) folgt nach Integration 
D) 
; ys ; , & 1+ : ' ‘ 
er) oo r(J xy dt) dr + Cy or + a ; We — a, dt+2C, 9+ Aro; (27) 
— 0 J 
: 6 


so daB fiir die Spannungen des Grundzustandes nach (9) 


E |d : 
C50. 4 Ler ae c—_ 5 (V Uo) — sae fo dt 
E [ (6 Ey (1 Py) (Lee y) I af r(f 28; 
= Jy) Cio +A, 9 2 Sate 2 | fie. dt) dr| , (28 a) 
E ity 
Go,0 = ips ~ - = oe AG Uy) == ao Oy, w 
0 
E 4. 4 Gan l—2») +) 0—2%) peel fee an)ar\|. 
ocr NS Got ¥ Aro 2 oe EEF x pay ole ? 
(28b) 
- 
oe | ns Ay 4 ey Seay Up) — aes = ac a, dt} = A,,E—E i a, dt + (6,9 + 5,9)  (28¢) 
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folgt. Aus den Randbedingungen o, = — p, fiir r = 1; und 6, = — Pp, fir T= iT, (wir wollen hier 
noch die Wirkung eines Innen- bzw. AuBendruckes p; baw. p, beriicksichtigen) erhalt man 


1 hen) 1 ae C sa0i(L— 27) pt Pak. 
C9 + 9.Ay9 =o OE) [ {fovea rt Faep Ley Be 3202) 


a Ue 


ui 


Te re rl ¢ /?8 r=Tq rita | Ps Pa 
lite ar | | r{ J acat) ar “+(14 ) ae alee ees | (29b) 


Br 2 ee hee 
il vr eh Pate 


ica) 
be 
i—) 


und damit ergeben sich schlieBlich die Radial- und Ringspannungen 


C.=— eu) <p (Gi 1 fr (fo i ar—|l meat ic (soy it | 


EO) p,ri |f{ra\ E(9) pyre iy 
E(8;) 2-73 ( | 1 E(0.) 72 — 3 [ i) | (30a) 
— — EO) peeiae ae tea 5 ade dros 1+(2) I fr ford ae 
OO Te U r2 — 7? r if ee r ae 
EO) pri ra\? | E(#) [Dn he Tae 
‘ E(;) r2—r? (7) . i ~ E(8_) 72 — r? I! Ms =) | j (30b) 


Zur Berechnung der Axialspannungen 0, nach (28c) benétigen wir noch eine Gleichung fir 
die Integrationskonstanten, die wir den Lagerungsbedingungen des Rohres an seinen Enden ent- 
nehmen. Wahrend sich im Falle des ebenen Verzerrungszustandes (e, = 0) 


Ay) =0 (31a) 


ergibt, hat man fiir den Fall, daB an den Rohrenden eine gegebene Axialkraft P angreift, 


NS, AE = 5270, ( GO, pdr bs 
(F) rj 


woraus mit (28c), (30a) und (30b) schlieBlich 


Pa eye (’ 
ae ar | r EI oh it) dry i) (Or,0 “+ Gp,0) dr 


4,9 a a = 
frEdr 

P rae is 
iat Tay [rE (fa at)dr 

= T; i / 2 1 Z o diva ie Pp; r? Oni 31b 
a ee : = [orl fo, dt\dr Le 

{°r Er poses: => ( OT BB) — By) ETP) 
hervorgeht. 


Fir die Berechnung des ersten Zusatz-Zustandes (k = 1) betrachten wir die Differential- 
gleichung (26b), aus der nach fortlaufender Integration schlieBlich 


= 1—2 1 = a 
rere EOL Le ary dr + Gyr + Sn, 


(32) 


Spas Ml = (l—2»)(1 + = - 
V Uy er (r v,) 4 Ay, =— hers 2 [fe dr +-2 Ch Ay 
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folgt. Die zugehérigen Spannungen erhalt man hiermit nach (9) zu 


= E \du v a 
aa = 745 | ag i) 
E(Q—2y»)[1—» [ Weale Bei Ca tayed G 
Te = far—i| r (fa dr) dr| Tt ph ae (33a) 
A E [o i ee 
ee |” aa lV i) 
ee aa ee is ae ge Ee Caen Ome tee, 
apes Were s eS 2 I fo dr + — 3 | r (J fo dr) ar a7 a, (ee si 2), (33 b) 
64 Al, b+ @,, +3,,); (33¢ 
/ 


Fiir die Integrationskonstanten gewinnt man aus den (nunmehr homogenen!) Radialspannungs- 
Randbedingungen o, , = 0 fiir r =r; und r = r, zunachst 


Cyitvydy,  (+r”Q—2») 1 


aay [2 Pfoael, 2 — fr fear) “4 > (Bta) 


eyo es ees 
i 21 i} y TF 


g- e 2 
Li, es lis 


= (=) 272 ji Ta l r=Tq 
C41 = wees »: 2 af pane te dr — | [Ufo dr) ar : | (34b) 


so daf§ damit fiir die Radial- und Ringspannungen der ersten Zusatzlésung 


2 Badal —(29] fea) fe 


AEE) Uae) pe pase om 


vj 


som ee EY] foarte) a] fa 
Salle! 1 i fr(ftoa ars fee (2) [or Liar) de tart Fan 


gefolgert werden kann. Zur Berechnung von G,, ist nunmehr noch A, , zu berechnen, wofiir im 
Falle des ebenen Verzerrungszustandes 


4,,=9 (36a) 


und im Falle einer vorgegebenen Axialbelastung P des Rohres an seinen Enden die (ebenfalls 
homogene!) Randbedingung 


oe dF = 2x fro,,dr=0 
gilt, aus der 
fr Gr + Gp,1) ar 


ee sae! Las oe (36b) 


hervorgeht. 
In entsprechender Weise werden nun die iibrigen Naherungen gebildet, so das also allgemein 


fiir die Spannungen des k-ten Zusatz-Zustandes mit 


= Vir) 
isos ae r,k—1 (37) 


Spannungen der Grund- 
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schlieBlich 
= HOGG a ee ee ee 
Ok = 78 78 cet, [fe-rar Sa ; 1 
1—2y/[{/ra\? 
« fae nell 
7, =720™m v r 
f [J fe .dr ar—|i os ) 
kg/em _ ; 
a < [ r( J Je-rdr)d ) (38) 
= 0 F 2 
a= dela) far ale 
1—2yv/[[/ra\? 
x J fords =a ey 
x Jr (J fi dr) Jar+ [1+ (2 )| 
x ( r( J fea dr) dr — (6 Cama A) | va) P) 
(38 b) 
O. 5 =s oe E =a (6, +6, T, a) (38 c¢) 
und 
Ai eatOs (39a) 
- 1 (Gr,k + Go, k) dr 
baw. A,, = —v ——___— (39 b) 
JrEd 
folet. Der endgiiltige Spannungszustand wird damit 
im Sinne von (9) 
GO, =G,.9 + De 6,5 0,0 1 Oa | 
Abb. 2. Wiarmespannungen. i 7 i ! 


lésung (k = 0) 


— | 1 ok GF +2. 
+— Spannungen bei unver- Oy i 0,0 ap » € Ck On.0 aie 5 Og, k? (40) 
dinderlichen Stoffwerten 1 
Ate %40)5 d n n 
pannungen der ersten val Ibo ~_ a 
Zusatzlésung (k= 1). = O20 at » € OL ma O20 i 2, O,, k 
(100fach vergroBert.) L il 


In Tabelle 1 sind einige Zahlenwerte fiir die Spannungen der Grundlésung k = 0 und der ersten 
Zusatzlésung eingetragen’. Man erkennt, daB die Funktionenreihen (40) so stark konvergieren, da 


12 
13 
14 
15 
16 
17 
18 


stimmt. 


°r,0 oy, 1 | or %,0 | Fp, 1 | 0 ee Gz 

[kelem'] | [kg/em'] | {kg/em*] | [kg/em') | [kg/em!] | [kg/em] | [kgjem*] | [kglem!]  [kg/om] 
— — 

E 14,0 —1499,9 | —1672,6 —1,3 —]499,9 

—101,6 | —0,6 | — 91,12} — 985,1 —12,8 | — 886,4 ] —1096,9 —I,1 — 976.5 
—144,6 —0,4 —128,77 | — 404,1; —11,5 | — 365,0] — 558,9 —0,5 — 493,8 
—145,6 —0,07 | —129,20 | + 102,0| —10,4 + 85,8] — 54,0; -+0,01 |— 43,4 
—116,7 +0,3 —103,18 | +. 548,7 | — 9,9 + 480,7 + 421,4|) +0,5 + 377,2 
66,3 | -+0,2 |— 58,32] + 947,2| —8.6 | + 831.4] + 870.5| L077 | 4+ 772-6 

0 0 | 0 +1300,3 | — 8,1 +1146,1 | +1295,8 | 0,9 1146,3 


,2 Der Zusatz-Zustand k = - 1 wurde hierbei durch numerische Integration der Gln. (35a, b) und (36 b) be- 
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man bereits die Grundlésung k = 0 als hinreichend genaue Naherungslésung ansehen kann. Die 
jeweils dritte Spalte in der Spannungstabelle gibt die entsprechenden Warmespannungen G6 nach 
der Theorie unveranderlicher Materialwerte an, die durchweg etwas kleiner als die genaueren 
Spannungen sind (Abb. 2). Die wirklichen Warmespannungen in einem Stahlrohr sind also gréBer 
als die nach der Theorie unveranderlicher Stoffwerte E, und oj; berechneten. Die Unterschiede 


B; =350°C B, =350°C 


7009 


400 


Abb. 3. Spannungen infolge rein mechanischer Einwirkungen. a) Belastung durch Innendruck p;; b) Belastung durch AuBendruck pq; 
Spannungen der Grundlésung (k = 0); —.—.— Spannungen bei unveranderlichen Stoffwerten. 


sind im vorliegenden Falle noch nicht allzu groB+ (rund 10 bis 15%); doch kénnen sie erheblich 
sein, wenn fiir das Temperaturfeld eine andere Randwertaufgabe vorliegt, wie vom Verfasser friiher 
gezeigt wurde. Die aus rein mechanischen EKinwirkungen herriihrenden Spannungen andern sich 


demgegeniiber nicht in demselben MaBe (Abb. 3). 


b) Zweites Beispiel. Als weiteres Beispiel untersuchen wir einen Stab konstanten Querschnittes 
(Abb. 4), an dessen Enden ein in der Symmetrieebene (x, z) des Stabes wirkendes Moment M) und 


Abb. 4. 


eine in derselben Ebene wirkende Langskraft S angreifen sollen. Der Stab stehe auBerdem unter 
dem EinfluB eines allein mit z veranderlichen Temperaturfeldes ) = /(z). 


1 Ermittelt man mit Hilfe der Theorie unveranderlicher, fiir die mittlere Rohrtemperatur %,, = 315° C 
giiltiger Stoffwerte die Warmespannungen, so ergeben sich hierfiir wegen E,, = E (315° C) = 0,891651 E, und 
Otm = 4 (315° C) = 1,2625 a0 [mach (22)] die Spannungen o’ = (Em Otm/E a0) @ = 1,125709 6, und diese 
stimmen mit den Spannungen der Grundlosung sehr gut tberein. Dies liegt daran, da die Temperaturunter- 
schiede innerhalb des Systems nicht sehr grof sind (max A? = G; — ba = 70°), so dab sich die tatsachlichen 
Stoffwerte nicht sehr von den mittleren Werten E,, und o;m unterscheiden. 
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Ebenso wie in der gewohnlichen Balkentheorie kommen wir auch hier zu einer Lésung der Ver- 
schiebungsgleichungen, wenn wir voraussetzen, dafs die Spannungen o, bzw. o, gegeniiber den (von 
y unabhangig angenommenen) Spannungen o, vernachlassigt werden kénnen. Sehen wir auBerdem 
auch noch samtliche Schubspannungen als vernachlassigbar klein an, so entsprechen diese hin- 
sichtlich der Spannungen getroffenen Annahmen auch im vorliegenden Falle, wie in der gewéhn- 
lichen Balkentheorie der Bernoullischen Hypothese?. 

Mit o, <a, sowie 6, <o, nimmt das Hookesche Gesetz in guter Naherung die Form 
a2) A(2) 
=> Ae} pat &, =E ag i 3 x, dt (41) 


6 ) 


E 


an, und da aus der Gleichgewichtsbedingung am Kérperelement in Axialrichtung 0¢,/0x = 0, also 
6, = 06,(z) hervorgeht, sind demnach auch die Dehnungen ¢,, ¢, und ¢, allein Funktionen von z, 
so daB mit den zunachst unbekannten Funktionen g,(z), go(z) und g,(x) fiir die in die (x, z)-Ebene 
fallenden Verschiebungen 


- = U, Cy aa u, €. = [x &,(2) a g,(2)| C., oe [go(z) te 83(x)| e. (42) 
folgt. Mit 
0 0 
BAS ea OF 
erhalt man hieraus 
= dgy = alee (ables oad Gt Gk ale, 
TN z)+t+—; Vi ae —- + —=|; VKuU=—el«*«— — SS 
7 U = e,(2) +33 / (Vu) a de) VX U ex BS z dx) 
as are Geis URS ‘d*g, dex 
AN] Wee (x dem ae} ARS es ze) 
be de, AB deen dg, 
Kibe== VV Wie (YX) (x aa a) + e, eS + de) 
und damit werden schlieBlich die Verschiebungsgleichungen (16) 
Peng , dg, 
ee 7 FEE ols (43 a) 
b 
d*g. 1—~» d’g,. 1 + v dex, al ff 
Eos qa ee?) aoe (43 b) 
i 
A d* Ex,1 d’g,, 
tee af =='04 (44a) 
ago. 1—v @g,,, 1 + v déx,y 
dz? 1 2 dx? 2 diaa le be (44) 
aéxk , CEik 4 
22 a5 dz — 0 r) (45 a) 
dgou | Lv deo le aoe . 
A eben ere el rg ae Gol 


Wir betrachten zunachst den Grundzustand (k = 0), also die Gleichungen (43a) und (43b). 
Gleichung (43a) wird befriedigt durch 


d’g,, . 

ae = 0 5 also £1,0(2) = Ajo z+ Ay 9 (46 a) 
sowle 

dex 

“a 9 also eg 0l2) = Cot + Cao, (46b) 


so daB damit Gleichung (43b) die Form 


6 
g9____1—vd¥,q 1+ 7,, avi 
d2 Tee eG) Cio + (lL + 7) a a, dt 


0 


* Wir haben es hier also mit einem speziellen [in der (x, z)-Ebene wirkend 
> 2)- en| eb S 
zu tun, so da wir nachfolgend die Verschiebungsgleichungen (14) baw. (16) iiauteeionee igi 
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annimmt. Fortlaufende Integration liefert 


dgy, = Gin. il o 
ee ( ae 7 G1o)# + (+9) fa dt + Cao, (47 a) 
1—» da 1 2 Cae 
s20(2) = —( las SAE eG +9) | (foudt de Ae Cae Cy ox GTB) 


Zur Reduktion der Integrationskonstanten betrachten wir die Hookeschen Gleichungen (41). Hier- 
nach muB 


o E o 
Bo ey eae 
\ é i} 


sein, und hieraus folgt nach Einsetzen von (46b) und (47a) mit einem Koeffizientenvergleich hin- 
sichtlich z* (k = 0, 1) 


. C30 =e Cy 0 
sowie 
dg, 
de + Go = 9, 
also 
x 
83,0(*) == Coz 5 Ag 9 x + Ay 9 ’ 
wobei A, . = — A, q ist, was sofort aus T,, 9 = 0, also 
dUx,0 Ouz,o 
=a oe ! 
baw. 


_ dex, 9 \ dg, | dg, 
a dz | dz | dx 0 


gefolgert werden kann. Damit lassen sich nun endgiiltig die Verschiebungen und Spannungen des 
Grundzustandes (k = 0) angeben: 


U, o(%X, 2) = ks (Cho z+ C, 0) “6 Ajo So Ag 4 ’ 


er 7d me \ 2 
tt, o(%> 2) = (1 +) | (J ay ds—9 (Cog fe Ca,92) — Chg dp nt Bys: (Bp = Ag gta) 
(48h) 


(48 a) 


/ F \ / 


oO \ 
0,,9(2) =E (éx0 — fo at) = E (Cra z+ Co,o— J x, at) 


Die verbliebenen Integrationskonstanten kénnen nunmehr mit Hilfe der Randbedingungen er- 
mittelt werden. Zunachst einmal folgt aus den Gleichgewichtsbedingungen 


ord 757 \o,2dF =M, 
(F) i 


(48 c) 


(F) 
mit den Abkiirzungen 
| de 
; 0 
tums =Usm = Fh oan es a ‘. 7 aa Pe (49 a) 
| nF | a ( | E* | 
(F) (F) (F) 
| zi 
ue = ———— =f2dF 49b 
tum = J, Ee res re ae , ik eee 5 ( ) 
| on ( Ey 
(F) (F) (F) 
E 2 
lene dF 
+ (F) 1 
oF Twp Eaar[[Enar\ aie 
(F) 
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. oO a oO ° 
E Er | ee | E | E oy 
Sz — —dt|\dk — | — 2d |) ——z — dt|\dF 
| Ey oy | all mo Ey J <£, Xo 
19) on Se eer Bs! a 


ae diag) (Oh ae yes J . . (49a) 
Stat E 
— 2 dF - z dF 
7 | gee (| im 
(F) (F) (F) 
a 
E E ot E f E ot 
= —— — dF — | —zdF =f ole 
| Ey Ps Ey | a J £, J Eo J) %€o 
apne Es) (Ns UE) Oo a oe (49 e) 
Ye a E , E 2 
Gl — 2° al zdF 
[a | ( pee) 
(F) (F) (F) 
fiir die Konstanten C, 9 und Cy 4 
S My Ay 1B 
Cio i: E, Fh *ms + E, Jy “MM = Fe ag i (50a) 
: Ss M. 
Coo = E, FSS — Ee PRYSM + O19 To xsr > (50b) 


so da® fiir die Spannungen der Grundlésung endgiiltig 
E S Myh J z 
Ox,0 = EF (xs *M S 7 Jy Per S Ey Ge i) 
4 ie N 


= o 
T2E 016 Ty (tgp | Es i 
010 Th h ; 


gefunden wird. Die iibrigen Konstanten folgen aus den Lagerungsbedingungen des Stabes. Setzen 
wir z. B. willkiirlich u, = 0 fiir x = 0 und z = 0 sowie 


x [u,, o(x, 2 = 0)] = 0 fiir Rea) 


d. h. an der Stelle « = 0 die Verschiebung und die Tangentenneigung der Achse z = 0 des Stabes 
gleich Null, was einer Kinspannung des Stabes entspricht, so ergeben sich 


Oo 
0 ee eee) | (fede) ae 
x 0 2=—0 


so dab der Verschiebungszustand der Grundlésung letztlich 


uz 9 —* (Cyo B= Cy, 0) 2 (52a) 
. o BP 2 
eee ») | [Jove dso (G05 ue C95) — Cs (52b) 
z—0 , 


wird. Insbesondere folgt fiir die Verschiebungen der Stabachse z = 0 


U,o(%, 2 = 0) = u(x) = Cy gx, (53 a) 


xe 


U, o(X, % = 0) = w(x) = — C, y ou (53 b) 


Wir betrachten nun die Zusatzzustande k = 1,2,...n, deren Differentialgleichungen (44) und (45) 
P 

homogen sind, also keine Stérfunktion enthalten [wiez. B. (43 b), wo der Ausdruck (1 +) d (jc) \dz 
0 | 


auftritt]. Hier werden also die Lésungen allein von den homogenen Lésungen der Differen- 
tialgleichungen gebildet, und da die zugehérigen Randwertprobleme simtlich homogen sein sollen, 
so sind also samtliche Zusatz-Zustande Null! Der Grundzustand stellt in diesem Falle die exakte 
Lésung des Problems unter den getroffenen Naherungsannahmen dar. Man erkennt nun aus (52a) 
und (52b), daB auch im vorliegenden Falle des temperaturabhangigen, also mit der Héhe z ver- 
anderlichen Elastizitatsmoduls die Bernoullische Hypothese vom Ebenbleiben der Querschnitte 
Giiltigkeit hat. Die Neigung du, /dz eines Stabquerschnittes x = konst. ist ebenso groB wie die 
Tangentenneigung dw/dx der verformten Stabachse z = 0, d.h. die Querschnitte stehen nach der 
Verformung weiterhin senkrecht auf der verformten Stabachse, und nach zweimaliger Differentiation 
von (53b) erhalt man mit (50a) 


d?w ; S M | 
= wy" =—C, 4 oo Ty 


ee ae Bee ee ix 
dx? EE Eh *Ms Ey Jy MM h XMT » (54) 
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d. h. eine Beziehung, die fiir konstanten Elastizitatsmodul, konstante Warmeausdehnungszahl und 
eine in der Balkentheorie iiblicherweise linear mit z veranderlich vorausgesetzte Temperaturverteilung 


H(z) = Ty + Ts 


wegen der in diesem Falle sich ergebenden Werte 


Fhag 2. F ie 
ire he = . tae Ninn = IL = tgp =Xyr = | 
MS SM ie Ss MM a Bee ST MT ° 
s 1 j 
Jgy = J, — BF, ts = zdF 
* . . (F) 
in die Differentialgleichung der elastischen Linie des Stabes 
i — My —S “Sm mel, Xto T; x d Msy to 1B 
E, Jsy h Ey Jsy h 
bzw. wenn die Stabachse z = 0 die Schwerachse ist, in 
aw —— My — a0" d 1 
EJ h 


tbergeht. Die Differentialgleichung der Axialverschiebungen erhalt man nach einmaliger Diffe- 
rehtiation von (53a). Es folgt mit (50b) 


du S My 
cae = Coo = ir “ss E, Fh Ysm + O49 ToXsr > (55) 
: : : eas 6 
und hieraus geht fiir FE = konst., «, = konst. sowie #(z) = T, + T, 7 die Relation 
a , My—S2zg a5 Msy ats) 1, %9 Ty T 
UI Hert 8 oBaje, O19 Ly E,F E, Jsy'S + 9 To ae aN ae (w F aa )ss + O19 £6 
S z 
=F pt sw" + 9 Ds, ( ge DT ysk Tt] 
und, sofern die Stabachse z = 0 die Schwerachse ist, 
S 
uw = E, F + O19 To 


hervor. Es erscheint daher gerechtfertigt, die Gleichungen (54) und (55) auch fiir veranderliche 
Schnittlasten M, = M,(x) sowie S = S(x) als giiltig anzusehen, womit man nunmehr das Problem 
der Stabbiegung auch im Bereich héherer Temperaturen allgemein erfaBt haben diirfte. 

Eine gewisse Vereinfachung ergibt sich in den Gleichungen (54) und (55) noch dann, wenn, 
ebenso wie bei den Verbundquerschnitten im Stahlbau und Stahlbetonbau, die Stabachse z = 0 
durch 
dE 0 (56a) 
Ey 

(F) 
definiert wird, womit xj,, = %s;4 = 0 verbunden ist, so daf die fiir diese Stabachse giiltigen Glei- 
chungen die Form 
a My X10 405 


fs cae aloe Ede RM Mem — po oN Te (56b) 
ee BE, F%SS + Oo Toxsr (56c) 

mit ; 

eer ee free NS ia 
‘MM ~ Fr 7 Ss iz 
— 2 dF = dF 
| 574 |, 
(F) (F) 


annehmen. 
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Bevor wir auf die Berechnung der x-Werte fiir ein konkretes Beispiel eingehen, sei noch ange- 
merkt, daB sich die durchgefiihrten Uberlegungen leicht auch auf den Fall der Doppelbiegung er- 
weitern lassen, sofern man fiir die Spannungen an Stelle von (48¢) 


’ oO 
0% = B (Gra z -|- Cy 0 L C3.0¥ an X, a 


einfiihrt. Die Ausdriicke werden entsprechend verwickelter, und man erkennt, daB auch die Haupt- 
tragheitsrichtungen von der Verteilung eines gleichzeitig vorhandenen ungleichmaBigen Tempera- 
turfeldes abhangig werden. Auf eine nahere Untersuchung dieser Erscheinungen wird jedoch hier 
nicht eingegangen. 

Eine weitere Bemerkung ist im Hinblick auf (51) erforderlich: Betrachten wir speziell 
den Fall eines statisch bestimmt gelagerten Tragers unter alleiniger Temperatureinwirkung, wo 
bekanntlich keine Schnittlasten M(x) baw. S(x) entstehen kénnen, so erhalten wir aus (51) 


ah Zz Ii 4 
6,.7(2) = E(2) % 9 To{ #7 + 7, “MT, fF, — dt (58) 
0 


also trotz des Fehlens jeglicher Schnittlasten Spannungen, die nur fiir x, = «,) und einen linearen 
Temperaturverlauf 


een 2 


verschwinden. In diesem Falle werden x 7 = %y7 = 1 sowie 


so da dann in der Tat o,,7 = 0 folgt. Wahrend die iibliche Rechnung mit unveranderlichen Stoff- 
werten also iiberhaupt keine Spannungen liefert, entstehen somit auch bei statisch bestimmt 
gelagerten Systemen Temperatur-Restspannungen. Sie 
bilden in jedem Querschnitt ein Gleichgewichtssystem, 
sind jedoch, wie das nachfolgende Beispiel zeigt, i. allg. 
vernachlassigbar klein (Abb. 7). 


An einem freien Tragerende klingen diese ein Gleich- 
gewichtssystem bildenden Restspannungen innerhalb eines sehr 
kleinen Randbereiches auf Null ab, wodurch sogar das Ent- 
stehen von Temperaturschubspannungen in diesem Randbereich 
bedingt ist. Diese sind jedoch wesentlich kleiner als die Tem- 

AbBYS: peratur-Restspannungen in einiger Entfernung vom Rande, 
so da’ der Spannungszustand im Randbereich i. allg. nicht 
von Bedeutung ist. Man erkennt diesen Sachverhalt am einfachsten an einem Trager mit schmalem Recht- 
eckquerschnitt (Abb. 5), in dem wir in erster Naherung einen ebenen Spannungszustand voraussetzen kénnen, 
und an dessen freiem Ende « = 0 wir den Randspannungszustand ox = — ox T(z) zusatzlich in Ansatz zu bringen 
haben, um insgesamt der Randbedingung der Spannungsfreiheit des Querschnittes x = 0 geniigen zu kénnen. 
Von diesem Randspannungszustand miissen weiterhin folgende Randbedingungen erfillt werden: 
h 


6, = 0 ’ te% = 0 ie -|- ae 


sowie 
Rc x 
==), C2) tiith ges @ ioe | > O., 
1 
Ein solcher, diese Randbedingungen sowie die Gleichgewichtsbedingungen des ebenen Spannungszustandes 
00z 


+ SF=0 


erfiillender Zusatz-Spannungszustand at sich mit dem zunadchst unbestimmten Faktor @ durch 
x 
Byes 


Ox(X, z) = —(1 +0 i) e * oy t(z) , (59a) 


afh 
x 7 ’ x 


‘ Ox T a :] =—o-o ; om 4b WHE) o (59b) 


2/h = —1/2 


9 == 
Tx2(X, 2%) = — Ww =e 


=e 
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a/h / 2[h 
6(x, z) = + ma —o 7 ee OxT a) d (7) 
s/h =—12 \2/h=—1/2 , 
a'h 2{h 
ee 2/1 C oe Ons g oxTd Shea : ae d sc 
1 @ 0) h e h x h i xT h 
2/h =—1/2 2/h—=—1/2 
BA ee 

= + o#(1 —o *) e BAT (2) (59 c) 


charakterisieren', und dieser kann nunmehr aus der Minimalforderung fiir die Formanderungsenergie er- 
mittelt werden, womit das Problem des Randbereiches naherungsweise gelést ist. 

Es ergibt sich schlieBlich, wenn man, wie in der Balkentheorie iiblich, //h S 1 voraussetzt, in guter Na- 
herung mit den Abkiirzungen 


4-1/2 ta or 4-1/2 Bei? 
a z z lla /ize OxT Y, 
Toke <T — mae ie) he asd Eat ae Nig “T £9 Zz 
es a ah | ea) em | es). 
== 1/2 Sail 2 Eesith 2 
die Formanderungsenergie 
l + h/2 
b iia 
W=> x lee + 2 — 20 ox oz + 21 + ¥) te] dx de 
x=0 s=—h/2 
eee [kp +2[1 +k ky] w? + hey wt 
ae Wi wae (1 + v) ka + v keg] w? + leg w*], 


und damit folgt aus dW/dm = 0 die biquadratische Gleichung 
2 k k. 
eee —_— - p coil — Ee 2 SS ee 
Coe G + v) ee y ale ATE 0 


mit der (einzigen positiven reellen) Losung 


i{ ieee ee Keo/keg 
o=|/za+opreg lV ee peer ik ee 


Auf eine genauere zahlenmaBige Berechnung dieses Abklingungsfaktors sei hier verzichtet, und es wird 
lediglich eine Abschatzung unter Zugrundelegung eines konstanten mittleren Elastizitatsmoduls E,, durch- 
gefiihrt. Setzt man fiir oz 7 die ein Gleichgewichtssystem darstellende quadratische Parabel 


OrT = 6) [1 — 1227], (¢ = 2/h) 


(Abb. 5) an, so wird 


P= oy (C— 40%), a dl gf? + 16¢4) = — 791 — 402° 
und 
_ % 4 pase ees ie aoe 
ty aE eerie porast 2 e640. 


so daB aus (61) (die Querkontraktionszahl fallt heraus) schlieBlich w = 5,03 hervorgeht. Damit sind die Span- 
nungen im Sinne von (59) berechenbar. Eine Extremalrechnung liefert an den Stellen 


es eae Le aie Sa gee 
P= a 0,2, Se er /3 = + 0,289 
fiir die extremalen Schubspannungen Ze 
o V3 
|Trzextr| = % oh w 0,356 a 
und bei 
Sega Zoe 
ow ee h 4 


fiir die gréBten Normalspannungen 


Ozmax = 09 16 2 ~ 0,214 Oo- 


Von diesen Extremwerten fallen die Spannungen dann mit wachsenden x/h rasch auf Null ab, was man schon 
daran erkennt, da® e—“+/h fiir x/h = 0,92 bereits nur noch 0,01 ist . 


Ee . h Fale : 
1 Diese Relationen liefern in .der Tat o, =Tx;=0firz = + a denn es sind Js oxT dz = 0 sowie 
—Nhje 


a hj2 . . - . oI . . 
J sox7 dz = 0, da die Spannungen ox 7 nach Voraussetzung in einem Flachenquerschnitt ein Gleichgewichts- 
—hjh 
system bilden! 

29 
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Wir kommen nun zur Berechnung der Beiwerte x nach (49), und zwar unter der Voraussetzung, 
daB die Balkenachse z = 0 gleichzeitig die Schwerlinie sei ( fed = . , fiir ein linear anwachsen- 
(F) 
1 


des Temperaturfeld! A 
== h eas Beet 62 
+e | ( 7) (62) 


sowie fiir die Materialgesetze | 


B(0) = By — By =F |1 — at (1+ Bz) =2,|1—<, (1+ Pz) (é = ol (63) 


1+ “27/1 + 65) 1 +26(1 fs Bx) | (6 = Tae} (64) 


Ato 


cA 
& 


z T. 
Hs) = T+ T,5 =T ! +Re 


O(9) = Xo + O41, F = %o 


= %9 


Es ergibt sich mit den Querschnittswerten 


wt i i ° Ie 
J, =f 2dF, Jaa | ear, Symae | AaF, ie =p | *4F, i, =V2 @) Fi 
(F 


) h 
F 
= (F) (F) ) 


sowie den Abkiirzungen 

TT’ HIN , . 2 i" * ye 
By =1—4,(14 2B + RP, B,=1—4)1+6(F) |, B= O87 (2 P5,)| 
2 Sy, s i iy ay 
JRO +46) —C.Sa(F) 7 | (66) 


Dy = (La) (L4 La) + BF] Co 80. 664 0) EB 
fe 


2 an 8 fm s es ; es e is 
he ee ISIS el MG I OT (1+4¢,) eee 
schlieBlich 


h 
Poe: Emig ie hie Ae ibe k 
G55 2 BiB, Bee | VSM em payee Be UM Dibea Bae (67) 
i h 
B,D, Sep B; D, a pita a 
Sie = B,B, = Be a B, B, — B “J 


Zahlenbeispiel: Wir betrachten einen Stahltrager IP 20, an dessen Unter- bzw. Oberkante die 
Temperaturen #,, = 350°C baw. 3}, = 250° Cherrschen sollen. Demnach sind T, = +(3, +0.) = 300°C, 
T, =0,— 0, = 100° C, und fiir die Stoffwerte hat man 

Ey =2,1- 10% kg/em?, Eo == 2, kelem* “C2 (200° C = 9 = 400 1G), 
Opp ise LOR Get. XO LO ate Gae3 


Es ist also 


pee br eee _ ET) - 1 Ty x 
coal e: = Eee 8 a a 
und die Querschnittswerte betragen 
F a 82,7 cm2, J; a 5950 cm‘, J = ae — (0) 5 
yy 1 b he d Ne 
——— 4 —= —_ ao ease EES SRO) pee, — 4 
Jy he | # dk 80 ! i 5)( 2 i) == 1205 em, 


(F) 
(b = Flanschbreite, t = Flanschdicke, d = Stegdicke). 
Mit diesen Zahlen erhalt man schlieBlich 


ug =1,101943, gp = Hyg = — 0,072 818, Xp = 1,102 218, 
Pep =O ae c= Th OCR e, 


und man erkennt hieran, daB sich z. B. reine Biegedeformationen an diesem Trager um (#7 — 1) x 
100 = 10%, die aus einer ungleichmafigen Erwarmung hervorgehenden Temperaturdurchsenkungen 
sogar um (x77 — 1) 100 = 25% gegeniiber denjenigen, die mittels unveranderlicher Stoffwerte 
berechnet werden, vergrdBern. Hier kénnen also schon merkliche Unterschiede gegeniiber der 
Rechnung mit unveranderlichen Stoffwerten auftreten. 


1 Das genaue Temperaturfeld ist mit Hilfe von (21b) ermittelbar, jedoch kann man dieses, sofern die 


Randtemperaturen, d.h, T, und T, vorgegeben sind, bei Stahltragern mit ausgezeichneter Genauigkeit durch 
(62) annahern., ; 
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Aus (51) ergeben sich damit folgende Spannungen 


_a) infolge alleiniger Langskraftbelastung 


SF =|! (1 +8 Z) 


b) infolge alleiniger Momentenbelastung 


nage 22D (E+ By) 


c) Warmespannungen 


ss —_4*ms i , 


Opp = Eq ayy Ty 


iol bbs) 


In Abb. 6 sind die Spannungszuwiichse 
Myz S 
SAG. ~ ee __ Ox,My 92 Ac ener oe) = 
M/W M/W MyjW “h’ ‘S/F S/F Ss 
gegentiber den iiblichen Balkenspannungen dargestellt. Man erkennt die sehr geringen Abwei- 
chungen, und sieht auch am Beispiel des Stabes, daB die Veranderung des Spannungszustandes 
infolge einer rein mechanischen Belastung bei Vorhandensein eines Temperaturfeldes nur sehr 


Ox, My —_ 


AG AG 
4/7 » SF 
4 =250°C 0,01 0001 


8=350°C 
Abb. 6. Spannungszuwichse / o gegeniiber den iiblichen Balkenspannungen. 
a) Biegebeanspruchung; b) Normalkraftbeanspruchung. 


8, =250°C 


B= 350°C 
Abb. 7. Wirme-Restspannungen. 


geringfigig ist. In Abb. 7 sind die Warme-Restspannungen dargestellt. Auch diese sind im vor- 
liegenden Falle nicht erheblich (sie kénnen gréfer werden, wenn zwischen Ober- und Unterkante 
des Stabes gréRere Temperaturdifferenzen herrschen). 

Zum AbschluB untersuchen wir noch das Problem der Langskraft-Biegung eines beidseitig ge- 
lenkig gelagerten Stabes bei gleichzeitigem Vorhandensein eines ungleichmaBigen Temperatur- 
feldes. Ausgehend von (54) erhalten wir mit 

M,(x) = P w(x), S(«) = — P 
die Differentialgleichung 


P P ro T. 
w+ tay ee — (rs E,Fh a a 4 a2) 
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mit der allgemeinen Lésung 


. 2 5 7 ; : ‘i fe pa 
w(x) =—1(%) oe ! “MT Hoa TE + Ccoshx + Csiniz) Gal a= iru ey) 


“MM “VM 


Die Integrationskonstanten werden aus den Randbedingungen w(0) = 0 und w(l) = 0 zu 


» 
“ 


M4 
MS 
C= (| Ss 4 


“mr Ey %o TF | i 
Pe a > 
“VM 


“MM iF 


*MS 


sin/jl 


ump Hy %o Ty 4 1—cosdl 
“MM vid 


\*“MM 


bestimmt, so da® damit letztlich fiir die Durchbiegungen 


Fone 2 Mera Lily Otay Te LEN esis re 
w(x) = (7p) Ee pa MT a is jl 


“MM “MM 


und insbesondere fiir den Biegungspfeil in Balkenmitte 


“MT 


sin Al 


4 sin A (l — x) 1 
sin A | 


2 l (z & *M Ss 
= w(—|=h|— acre 
J (=) q) ee 


“MM 


1 ' 
x = 0,907 —_____—- 
Hy” 

400 ze 


Abb. 8. Lingskraft-Biegung cines ungleichmiBig erwirmten Stabes. 


genauerer Biegepfeil; 
—.—.— Biegepfeil bei unveriinderlichen Stoffwerten (M9, «;9)- 


Sar jt eB) 


iP 
cos 7: 
hervorgeht. Fiir P = 0 und unverander- 
liche Stoffwerte folgt daraus 
1 1\2 
fr, = 3 “t0 T; (i) h. 

In Abb. 8 ist 7 = f/fr, in Abhangigkeit 
von P/P),,, wobei Po,, = 2? E, J,/P die 
kritische Last fiir unverdnderliche Stoff- 
werte bedeutet, graphisch dargestellt (es 
wurden hierbei die Zahlen des vorangegan- 
genen Beispiels benutzt). Man erkennt, daB 
die Durchsenkungen mit wachsendem P 
schneller anwachsen (sie sind schon fiir P = 0 


um 100 (#77 — 1) = 25% gréBer als die in 


iblicher Weise berechneten), und daB sich 
die kritische Last verringert. Hier sind die Unterschiede zumindest bei den Verformungen und 
damit fiir die Momente und Spannungen schon kaum noch vernachlassigbar. 


5. Zusammenfassung. Die in Ziffer 4 untersuchten Beispiele zeigten, daB sich bei Vorhanden- 
sein eines ungleichmafigen Temperaturfeldes im Bereich héherer Temperaturen nicht nur die reinen 
Warmespannungsfelder gegeniiber denjenigen, die unter der Voraussetzung unveranderlicher Stoff- 
werte berechnet werden kénnen, dindern, sondern daB auch eine Veranderung der durch rein mecha- 
nische Einwirkungen hervorgerufenen Spannungsfelder auftritt. Diese lassen sich bei Systemen 
aus Stahl i. a. (sofern nicht extreme Temperaturunterschiede im System vorhanden sind) vernach- 
lassigen bis auf die Klasse derjenigen Probleme, bei denen die Spannungen am verformten System 
von denjenigen am unverformten System wesentlich verschieden sind (z. B. bei der Langskraft- 
Biegung). In diesen Fallen kénnen sich auf Grund der gegeniiber der Theorie unverinderlicher 
Materialwerte gréBeren Verformungen auch erhebliche Anderungen der Spannungszustande ein- 
stellen, so da’ man hier die Temperaturabhangigkeit der Stoffwerte beriicksichtigen sollte. Die 
Anderung der reinen Warmespannungsfelder gegeniiber denjenigen, die mittels unveranderlicher 
Materialwerte bestimmt werden kénnen, ist in der Regel, d. h. wenn die Temperaturunterschiede 
innerhalb des Systems nicht extrem grof sind, noch nicht allzu gro8, so daB man sich auch im 
Bereich héherer Temperaturen mit der Theorie unveranderlicher Materialwerte E und «, begniigen 
kénnen wird, vor allem dann, wenn man die einer mittleren Systemtemperatur entsprechenden 
Werte benutzt. Dies ist aber, wie der Verfasser schon friiher gezeigt hat, i. a. nur dann vertret- 
bar, wenn man der Spannungsberechnung das genaue, die Temperaturveranderlichkeit der Warme- 
leitzahl beriicksichtigende Temperaturfeld zugrunde legt, da dieses in manchen Fallen von dem mit 
konstanter Warmeleitzahl berechneten merklich abweichen kann. 


(Eingegangen am 4, Marz 1958.) 


Anschrift des Verfassers: Priv. Doz. Dr.-Ing. Rudolf Trostel, Techn. Universitat Berlin, Lehrst. f, Mechanik 
(Prof. Dr.Ing. I. Szabé), Berlin-Charlottenburg 2, Hardenbergstr. 34. 
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Die Kreiszylinderschale mit kreisférmigem Ausschnitt 
unter Schubbeanspruchung* 


Von D. Withum 


1. Kinleitung. Wird ein diinnwandiges Kreisrohr an seinen Enden durch Torsionsmomente 
belastet, so erleidet seine Wandung eine reine Schubbeanspruchung. Dieser gleichformige Span- 
nungszustand erfahrt in der Nachbarschaft eines Loches (Abb. 1) eine drtliche Stérung, die mit 
erheblichen Spannungsspitzen verbunden ist. Bereits 1898 wurde dieses Problem von G. Kirsch} 
fiir den Sonderfall der unendlich ausgedehnten Scheibe, also bei verschwindender Schalenkriim- 
mung gelést. Danach steigen am Lochrand die Hauptspannungen auf das Vierfache an. Zahlreiche 
Versuche® zeigten, daf} bei endlichem Schalenhalbmesser jedoch wesentlich héhere Spannungsspitzen 
auftreten kénnen, so daf} eine genauere theoretische Untersuchung wiinschenswert erscheint. 


4 


Ke 


Abb. 1. Kreiszylinderschale mit Ausschnitt Abb. 2. Koordinatensysteme der Kreiszylinderschale. 
unter Torsionsbelastung. 


Kine solche wurde von D. Thoma und M. Schilhansl? mit Hilfe einer Naherungsrechnung durch- 
gefihrt. Infolge zu grober Annahmen ergaben sich jedoch keine einwandfreien Ergebnisse. Wah- 
rend der Abfassung der vorliegenden Arbeit wurde ferner ein Aufsatz von A. I. Lurje* bekannt, der 
die auf Zug beanspruchte und durch ein Kreisloch geschwachte Zylinderschale behandelt. A. I. 
Lurje geht von der gleichen Differentialgleichung aus, die auch dieser Arbeit zugrunde liegt, wahlt 
jedoch einen anderen Liésungsweg und beschrankt sich auf sehr kleine Lécher. 

Im folgenden wird das Problem unter Beriicksichtigung der Biegetheorie der Kreiszylinder- 
schalen mit Hilfe der Stérungsrechnung gelést. Fiir sehr kleine Lochdurchmesser wird die Lésung 
in geschlossener Form dargestellt. Anhand einiger Zahlenbeispiele werden die Konvergenz des 
Verfahrens nachgewiesen und die Spannungszustande bei Schale und Scheibe miteinander ver- 
glichen. SchlieBlich wird noch ein Diagramm angegeben, das bei beliebigen Schalenabmessungen 
unmittelbar die maximale Spannung am Lochrand abzulesen gestattet. 


2. Rechnungsannahmen und Voraussetzungen. Es gelten die iiblichen Annahmen der Schalen- 
biegetheorie. Ferner wird vorausgesetzt, da die vom Lochrand ausgehenden Stérungen bis zum 


* Gekiirzte Fassung einer von der Fakultat fiir Bauwesen der Technischen Hochschule Hannover ge- 
nehmigten Dissertation. Referent: Prof. Dr.-Ing. habil. A. Pfliiger, Korreferent: Prof. Dr.-Ing. W. Zerna, — 
Der Verfasser ist den Referenten fiir wertvolle Anregungen und grofziigige Unterstiitzung zu Dank verpflichtet. 

1 G. Kirsch, Z. VDI 42 (1898), S. 797. 

2 4, Thum und O. Petri, Steifigkeit und Verformung von Kastenquerschnitten, VDI-Forschungsheft 409, 
Beilage zu Forschung auf dem Gebiet des Ingenieurwesens, Ausgabe B, 12 (1941); W. Bergmann, Grundl. 
d. Landtechnik, Heft 4, S. 12, Diisseldorf 1953; F'. R. Schlechte und R. Rosecrans, Experimental Stress Ana- 
lysis of Stiffened Cylinders with Cutouts — Pure Torsion, NACA TN 3039, Washington D.C. 1953; H.G. 
McComb, Stress Analysis of Circular Semimonocoque Cylinders with Cutouts, NACA, Washington D, C. 1955; 
P. Kuhn, Stresses in Aircraft and Shell Structures, New York Toronto London 1956. 

3 D. Thoma und M. Schilhansl, Luftf.-Forsch. 19 (1942), S. 210. 

4 A, I. Lurje, Prikladnaja Matematika i Mekhanika, Neue Serie 10, 3, Moskau 1946, S. 397. — Auszug 
in deutscher Sprache im Buche von G. N. Sawin, Spannungserhéhung am Rande von Lochern, S. 130, Berlin 
1956. 
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Schalenrand hin abgeklungen sind, und daf bei geschlossenen Rohren die Randstérungen sich nicht 
iiber den Zylinderumfang hinweg gegenseitig beeinflussen. 


3. Biegetheorie der Kreiszylinderschale. a) Geometrie. Zur Beschreibung eines beliebigen 
Punktes P (Abb. 2) der Schalenmittelflache wird ein orthogonales Koordinatensystem verwendet, 
das bei Abwicklung des Zylindermantels in die Ebene in ein Polarkoordinatensystem ibergeht. 
Mit den iiblichen Zylinderkoordinaten x, y ist dieses Koordinatensystem (0, p) durch die Trans- 
formationsformeln 

Xx == 0-COS Y, 

y =osnyp 
verkniipft. Dabei bezeichnet R den Schalenhalbmesser. Die Koordinatenlinien ~ = konst. sind 
Schraubenlinien der Schalenmittelflache. Der Lochrand wird durch die Koordinate @ =r be- 
schrieben, und der Lochmittelpunkt fallt mit dem Ursprung 0 zusammen. 

b) Differentialgleichung. Werden in der Schalenbiegetheorie gewisse Vereinfachungen vor- 
genommen, die u. a. von A. E. Green und W. Zerna begriindet wurden, so laBt sich das Problem 
auf eine einzige partielle Differentialgleichung vierter Ordnung fiir eine komplexe Spannungs- 
funktion p zuriickfiihren. Aus dieser kénnen nach Abschnitt c) die SchnittgréBen und die Radial- 
verschiebungen in einfacher Weise ermittelt werden. Nach W. Zerna? gilt 


| (ae R<psing S +28) (1) 


4 
2, Sl aD) 
AAy —2i02°% 0 mit i (2) 
: 
Hierbei bezeichnet ; - 
VY) pe eee 
Vis One dy? 
den Laplaceschen Operator, i die imaginare Einheit, 4 die Querkontraktionszahl und t die Wand- 
starke der Kreiszylinderschale. 
Wie sich zeigen laBt, ist y gegeniiber Koordinatentransformationen invariant. Wendet man 
die Transformation (1) auf die Gleichung (2) an, so folgt: 
: Oy 1 op 1 @y 2 ay OU \ ee 
2 gS Nee pas = 
SoA Ct (se Q@ 9% 7a) aa? 3 dg ap ee) sin2 | — ©) 
wobei der Laplacesche Operator die Form 


Anne Cpe oye? 
dg? ' @ AQ ' g? ay 
annimmt. Wird nun noch eine dimensionslose komplexe Veranderliche 
paeany Aiea oat 


eingefiihrt, und definiert man 


A*y= 


! oe 
so lautet die Differentialgleichung (3): 


A* (A*y 4 _ (Oy 1 op 1 ey\ ZO) 2 op\ . at 
(A*p + y) 4 (Fe E a ag cos 2 gy — E aap = ay) cindy =O. (5) 


% \ Mop 
Te 
TM, 
%p ? 
™, 
Mp Vo 
™M, 


™, 


} Abb. 3. Schalenelement mit SchnittgréBen. 
a) Liingskrafte und Schubkrafte, b) Momente und Querkrafte, 


’ ¢) SchnittgréBen und Formanderungen. Aus Abb. 3 sind die positiven Richtungen der 
chnittgréBen ersichtlich. Die Verschiebungskomponenten langs der positiven Koordinatenrich- 
tungen 0, y und z werden mit u, v und w bezeichnet. 


1 A, E.Green und W. Zerna, Theoretical Elasticit Sao i 
2 W. Zerna, Ing,-Arch, 20 (1952), $.357.. Pea 
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Setzt man mit dem Elastizitatsmodul EF 
LE Cyne Sea 


Ee pee die aus der Scheibentheorie bekannte Airysche Spannungsfunktion @ ein, so gilt nach 
. Zernat 


w=Rey und O= S Oily . (6) 


Hieraus ergeben sich die SchnittgréRen zu 


erect ee a (7a) 
ne Sor a , (7b) 
Py eall aun ae (76) 
m, = — BH (5s tee ge tn ar (7d) 
mya BE (WO Se tex + ae). (7) 
ee ay | | (7!) 
Mit 
B= pao 


ist dabei die Biegesteifigkeit der Schale gekennzeichnet. Da die Schubkraftverformungen vernach- 
lassigt werden, so sind die Querkrafte nur iiber die Gleichgewichtsbedingungen mit den Schnitt- 
momenten und damit auch mit der Radialverschiebung w verkniipft: 

" 0(Aw) a om 


Coe Bo je Tae ioe? (7g) 
a a (Aw) 1 « oM 
p= B E Op Le = en ~ Mee) 
Hierbei ist zur Abkiirzung 
m,—m, =—(1+ 4) BAw=M (8) 


gesetzt. 

Da mit der Differentialgleichung (5) an jedem Schalenrand nur vier Randbedingungen erfillt 
werden kénnen, werden in tiblicher Weise fiir den Rand € = konst. das Drillmoment und die Quer- 
kraft zur Randquerkraft 

x OMe 


Iors=16 E ap (9) 
zusammengenommen, 
Setzt man ferner 
so stellt 
i OG 
Q— ae 
2=2R(N4 ‘tM (11) 
eine weitere komplexe Spannungsfunktion dar, wobei entsprechend (8) die Abkiirzung 
N=n,+n, =A® (12) 


eingefiihrt ist. 
Die Verzerrungen der Schalenmittelflache sind durch das Hookesche Gesetz 


il 
ime Tar (Weal) Th.) | 
tf 
a a (ST (13) 
2(1 + w) 
6 “is op 


ror Siehe FuBnote 2 von Seite 436, 
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mit den Schnittkraften verkniipft. SchlieBlich bestehen zwischen den Verschiebungen der Schalen- 
mittelfliche und den Verzerrungen die Beziehungen 


ou Ww 


wo2R 


é (1 — cos 2), 
a 
(pt) tap ll +0829). (4) 


€ = & 


du Ov w . 
= add ; Ey se iy Sa y 
reo = E56 0+) 4 Rp sin2e: 


Damit haben wir alle erforderlichen Formeln zusammengestellt, um die Zylinderschale mit kreis- 
formigem Ausschnitt behandeln zu kénnen. 
4. Spannungs- und Formiinderungszustand der ungeschwachten Schale. Die bei der unge- 
schwachten Schale auftretenden ZustandsgréBen werden durch einen Querstrich gekennzeichnet. 
Das Drillmoment D (Abb. 1) ruft in der geschlossenen Schale einen konstanten SchubfluB 


D ° 

P= on Rt’ 4 

bezogen auf die Zylinderkoordinaten x, y (Abb. 2), hervor. Entsprechend verschwinden im Koordi- 
natensystem (0, y) simtliche SchnittgréBen mit Ausnahme von 


n =+Tsin29, | 
i, =—Tsin29, (16) 
Nog= + Tcos2@. 


Fir die Airysche Spannungsfunktion ergibt sich 


6 =— + sin2y, (17a) 
wahrend 
(e=T (17b) 
ist. Fir die Verschiebungen gilt 
u +H = Tsin2g, | 
é ( 
en tla as f 18 
mal + cos2q), (18) 
wi) 


5. Das Randstérungsproblem der geschwichten Schale. Ist in der Schale ein kreisférmiger 
Ausschnitt vorhanden, so kann der Spannungs- und Formanderungszustand aus zwei Anteilen 
zusammengesetzt werden. Der erste Anteil besteht aus den Beanspruchungen und Verformungen 
der ungeschwachten Schale, wahrend der zweite von einer Randbelastung des Loches hervorge- 
rufen wird. Diese besteht aus einer Gleichgewichtsgruppe, die den in (16) angegebenen Schnitt- 
groBen das Gleichgewicht halt. Die zugehérigen ZustandsgréBen seien mit ~ geschlangelt iiber- 


strichen, Die endgiiltigen Schnittgré8en und Verformungen kénnen dann durch Uberlagerung er- 
mittelt werden. 


a) Randbedingungen. Der Schalenausschnitt werde durch die Koordinate & = Er = Tabes 
grenzt. Dann gelten die Randbedingungen 


n, =—Tsn29, m, =0, 


‘ fi as 
ea ee fo ur 6 ==, (19) 


Fir omg co miussen die infolge der Randbelastung (19) auftretenden SchnittgréBen nach dem 
De Saint Venantschen Prinzip verschwinden. 

a Fourieranalyse. Die Zylinderschale ist symmetrisch in bezug auf die Koordinatenlinien 
gy =0 und » = 2/2. Der Spannungszustand ist fiir diese Koordinaten antisymmetrisch und mit 


periodisch. Die Langskrafte n, miissen daher in folgender Form in eine Fourier-Reihe ent- 
wickelbar sein: 


o 
n= Sin, sindng. (20a) 
n=1 
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Entsprechend lassen sich die ZustandsgréBen No» My, My, q,> Ia. N, M, &,, é,,U, w und die Span- 
nungsfunktionen darstellen. Fiir die Schubkrafte n,, ergibt sich eine in bezug auf ~ = 0 und 
gp = 7/2 symmetrische Entwicklung: 

ioe) 


Nog = D, Mpg, COSIND. (20b) 


Die ZustandsgréBen m,,,, dy> Yow kOnnen in gleicher Weise entwickelt werden. SchlieBlich gilt fur 
die Verschiebungskomponente v die Reihendarstellung 
s Lee a 
v= +%+ Dd’ v,cos2ngy. (20c) 
, itl 


Mit diesen Ansatzen wird die partielle Differentialgleichung (5) in ein gekoppeltes System ge- 
wohnlicher Differentialgleichungen vierter Ordnung 


s = Ll d*n—1 An—3 djr—1 , 2n(n—1) ~ 
* * 
 dintt  4n+3dyn41 2n(w+ 1) -~ (21) 
! 2 dt ae eet sdica . a = OU a ae 
j= 0 | 
umgeformt, wobei der Laplacesche Operator die Form 
A dp, 1 din An? . 


A*p, = Q,, == 


de ! & dle C2 Vn 
annimmt. 

Es ist zweckmabig, ein entsprechendes Gleichungssystem fiir die Funktionen Q,, aufzustellen. 
Wird (22a) nach #, aufgeldst, so ist 


py, =O fem fom dedl. (22b) 
Setzt man diesen Ausdruck in (21) ein, dann verschwinden die Doppelintegrale, und es folgt 
A*Q, + 2, +5 Qra1—n— lye fo4Q, dl | 
{ 
1 ee 
+ Quy tien +10 fo 410,146 =0 (n= 1,2,3,...),-— 9) 
Q,=0. 


Fiihrt man die Ansatze (20) in das Differentialgleichungssystem (14) ein, so erhalt man aus 
einem Koeffizientenvergleich 


: a l 1 l 
Cp OU oe (-, a ae tn] ; (24a) 
s ob 1 1 1 
Con vee (— 2 n Vv, a u,,) + 9 R («, | 9 Wn—1 | ®) Wr ) > (24b) 
y aad v : 24) 
Voon as ae (2 nu, —v, fia é v,) a IR (Way — W,—1) p) ( Cc 
wolyet Ti == 0h so 83 ona 
Eo) = En) = Vong = Ho = M = 9 
und 
Lee ame sic 


zu setzen ist. Ableitungen nach é sind durch einen Strich gekennzeichnet. 


Aus den Randbedingungen (19) wird 


Z == 9 tiie 71 3 ass fl Anite, == 1 | 
n = > 1 = r) 
f 0 fir n>1 ae 0 fir n>1f - fir€=&. (25) 
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Werden die SchnittgréBen durch die Spannungsfunktion @ und die Radialverschiebung w aus- 
gedriickt, so ergibt sich unter Beriicksichtigung der Reihenansatze (20) nach einigen kleineren 
Umformungen 

Phi ge IE Gite Tail 

| 0 fir n> If” 

Bape (2 tue) Es w, — [8 n+ 1—2np(2n—)] & Wy 


Eo P,, ap 2 nN P,, a 


So 5 ¢ 
+4n?(3 +4 (2n—1)] wv, =0, fies er 26) 
£ oo ile 7 = Il 
pte = te = ; 
() inkie 7SS dl 
£) we + hE Wn — fs 4 nw, = 0 


Fir €—> co muB jedes Reihenglied der SchnittgréBen einzeln verschwinden. 

c) Stérungsrechnung. Wird der Lochrand 9 =r einer unendlich ausgedehnten Scheibe 
durch eine mit sin2nq veranderliche Belastung beansprucht, so ist der Spannungszustand in 
der Nachbarschaft des Loches ebenfalls mit sin2 ng bzw. cos2nq veranderlich. Dies rihrt 
daher, daB bei der rotationssymmetrischen Scheibe keine Richtung g geometrisch bevorzugt ist. 
Bei der Zylinderschale geht die Rotationssymmetrie verloren, die Hauptkriimmungslinien sind 
ausgezeichnete Richtungen. Mathematisch driickt sich dies auch im Gleichungssystem (21) aus. 
Die Differentialgleichung vierter Ordnung fiir die mit sin 2 np veranderliche Spannungsfunktion yp, 
ist gekoppelt mit den ,,Nachbarfunktionen® »,_; und y,,1. Diese Kopplung ist um so geringer, 
je gréBer der Schalenhalbmesser R ist, im Sonderfall der Scheibe ist sie nicht mehr vorhanden. 
Fiir kleine Lochhalbmesser iiben daher die Anteile der Nachbarfunktionen keinen grundlegenden 
EinfluB auf die Gestalt der Lésung fiir ~, aus. Nach den Regeln der Stérungsrechnung werden 
deshalb diese Glieder mit einem Stérungsparameter ¢€ multipliziert : 


Aaa = iI ayn—1 4n—3 dijn—1 2n (n = 1) A 
* /A* 
A (A Vn Yn) é (3 dé? om (ee RE : dé é2 pil 
Ll @insi  4n+3 djn41 , 2n(n+1) -~ \ f (27) 
line 5) HE {= ; soe e | = oe ) Wn «| == (( (n = il oe ay oe 2) > 


4 


Wo == () 


Fiir ¢ = | geht (27) in das urspriingliche Differentialgleichungssystem iiber. Fiir ¢ = 0 entkoppelt 
sich das System, und die hierfiir zu ermittelnde erste Naherungslésung wird um so besser mit 
der wirklichen Lésung iibereinstimmen je kleiner |C| ist. 


Die Spannungsfunktion ~ wird nun in eine Potenzreihe nach ¢ entwickelt. Fiir jedes #, gilt 
dann 


i) n= 6h py) = (28) 
v=0 


Dabei soll pl) die Randhedingungen (26) erfiillen. Die p) fiir » > 0 diirfen also nur noch die 
zugehorigen homogenen Randbedingungen befriedigen. Auch die SchnittgréBen, Verzerrungen 
und Verschiebungen lassen sich in dieser Reihenform darstellen. 


Wird (28) in das Differentialgleichungssystem (27) eingesetzt und ein Koeffizientenvergleich 
fiir gleiche Potenzen von ¢ vorgenommen, dann ergibt sich 


2571) a Ge— > 
xy (v ~(y Li a Ay 8 aye Ss ) ue 
RAR oe ON es ee ee n—l 21 (ta-1) so 
A (A Gene aoes ) 9 dé? Ma me vile vex ¥ ay | 
PA 1) > Ane—D) 

ee Oe She Ty) (29) 

Tae ee ae a ee 
pe) = (0, pr <=). 


In dieser Differentialgleichung sind p—!) und pe) bekannte Funktionen. Die vollstandige 


Lésung fiir ~ setzt sich somit aus einer Partikularlésung 7) und der Lésung der homogenen 
np 


Differentialgleichung zusammen. Sie lat sich mittels Besselscher und Hankelscher Zylinder- 
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funktionen erster Art in geschlossener Form angeben. Werden mit a) , b) , ce), und d(), will- 
nes 5 4 = n, ti n, n, n, 
kiirliche Integrationskonstanten bezeichnet, so ist “ey 


yy (? vy) Fan v 1 v laa v fa i G 
PO) = a), C+ BO), + oy Jan(C) + a), HE () + GO. (30) 


nn €2n “ny 
Da die SchnittgréBen mit wachsendem Argument abklingen miissen, ist 


at?) = cr) =0. (31) 


n, 


Die beiden restlichen komplexen Integrationskonstanten gestatten, die vier Randbedingungen (26) 
zu erfiillen. 


Bei Kenntnis der jeweiligen Partikularlésung lassen sich die Funktionen pe aus der Differen- 
tialgleichung (29) sukzessiv berechnen. Da p = 0 ist, und die Randbedingungen (26) fiir n > 1 


homogen sind, ergibt sich fiir den Rechengang das folgende Schema, wobei die identisch ver- 
schwindenden Funktionen nicht aufgefiihrt sind. 


Der Ansatz (28) liefert in entsprechender Weise fiir das Gleichungssystem (23) 
1 ; - 5 
Be ey (2 n— 1) oe. fC 8 OC de | 
1 fon -—2n— v— f UG 32 
== MP — On $1) o fE-P AIO) dl (n=1,2,3,...), (¥=0,1,2,-.)s (ee) 
2O= 0, Q-Y=0. J 
Als Lésung dieser Differentialgleichung ergibt sich 
21) = 0, HIN) +90, (3) 
wobei nur der abklingende Lésungsanteil beriicksichtigt ist und 20) wieder eine Partikularlésung 


bezeichnet; s,,9 ist eine komplexe Integrationskonstante. 
Die Partikularlésungen der Gleichungen (30) und (33) lassen sich ebenfalls in geschlossener 
Form angeben. Das Ergebnis sei hier vorweggenommen. Es wird 


AC —— é v A = (ry) FP FW é 
I= TWO, de + Sd, LP HO (0, (34a) 
p=) p=0 
n l v F 
WS > pO HE ©) - (34b) 
p=2 p=9 
Die Randbedingungen fiir §— o sind erfiullt, da 
(ey on He, (ey) 0. (35) 
Zwischen den Koeffizienten der Gleichungen (34a) und (34b) bestehen die Beziehungen 
eee et Pep 1) 10), a) 


so), = — do), —2(p+1)d0,_,. (36b) 


n,p-+-l 
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Dabei ist 
abel 1s hi) aca d) 5 (”) On | 
0, p 0,p 0,p 0,p 
ho) =0 fur p= und p>7n, 
n,p Fi 7 (37) 
r?’) =—0 firs op ade and pe ae 
ae. =s’) =0 fir p< 0 id gpieny:s 
n,p n,p 


b@ und d©), sind als Integrationskonstanten frei verfiigbar. Die iibrigen Koeffizienten ergeben 
n,n n, 


. Gus = (y—I ei Moe 
sich aus den Funktionen p@—) und Mey! » baw. Q0— und 2°77 zu 


: nt pl ay) a ee Pe ey 8: 
be) =4(n—p+1)(n+p—1) Be) a + ane 2 be + Ee ON ip? (38a) 
. l pa (v—l ws i (v—1) 38b 
Ne es 4p (der) 4 oe ea) 2 Sn qe ( ) 


Die Formeln (36) folgen aus einem Koeffizientenvergleich, wenn (34) in die Funktionalbeziehung 
(22) eingesetzt wird. Die Gleichung (38a) laBt sich ebenfalls leicht beweisen, wenn (34a) in (29) 
eingesetzt und hier wiederum die Koeffizienten verglichen werden. Entsprechend erhalt man die 
Formel (38b) nach einiger Zwischenrechnung aus der Differentialgleichung (32), wenn dort die 
Lésung (34b) eingefiihrt wird. Damit ist die Richtigkeit der Ausdriicke (34) nachgewiesen. 

Die Koeffizienten dt) lassen sich auch explizit angeben. Aus (36b) und (38b) folgt die be- 


sonders fiir n = 1 zweckmafige Formel 


re ~ pat Pie Phe ayant 39 
dt) art ple 4p st F ee i > = 2) Pp! dae Lt” ( ) 
(aay 
Werden die Konstanten 
Pre ane (50) 5 = 30) 
be i is +- U Be p) dip a daa + U disp > (40) 
1’), = Te aE U rey 2 sv) == 3, ae ia 


in ihre Real- und Imaginarteile zerlegt, und setzt man zur Abkiirzung 
vie [ee ea =o Ze te) 2 Ot [2 Hy_,(0)] a Z,,p(é) ? (41) 
so folgt aus (34) mit (6) und (11) 


we) os ay (4 Ae cos eS se bt, sin a E2p Sle a, Zerle) > dv Zu pl) y) (42a) 


pP 


AUD) sa [085 e AG w\ 1 , a) ¢ 4) ‘ 
go —) DS [e by, sinh + 60), cos eS) ay + a, Zap(6) +d, Zns(8) ,  i(A2b) 
P 


») Et <(v We) ae aad! —() 7 __ a) & 
No = rR > (+ 7 cos sleet popieie a Hp tai Sn. p Lepp\e) = Sa ep oie (42c) 
P 
v 1 -b —(y . 7 aA(y x i —(y # Ay ps 
Mo) — 5 ey ( 7, sin i = al cos oe 2p aie ss Zn, p(E) + ae Linley (42d) 


Die Ausdriicke (41) ieee sich aus den Funktionswerten der Hankelschen Funktionen nullter und 
erster Ordnung vom Argument & yi unter Beriicksichtigung der Funktionalbeziehungen der 
Zylinderfunktionen berechnen. 

Die Randbedingungen (26) sind so gewahlt, da® sich die vier Gleichungen weitgehend ent- 
koppeln, wenn die Ausdriicke (42a) und (42b) eingesetzt werden. 

d) Konvergenz. Ein strenger Konvergenzbeweis der Stérungsrechnung diirfte nicht einfach sein. 
Im Rahmen dieser Arbeit wird daher die Konvergenz an Hand von Zahlenbeispielen nachgewiesen. 

Abb. 4 zeigt den Verlauf der Normalkraftsumme N() fiir verschiedene Naherungsschritte bei 


einem Lochhalbmesser & — 3,81. Wie zu erwarten war, ist die Konvergenz am Lochrand be- 
sonders gut und nimmt fiir wachsendes € ab. Ein ahnliches Verhalten wurde auch fiir andere 
SchnittgréBen nachgewiesen. 

Im allgemeinen wird es zu miihsam sein, den Verlauf der ZustandsgréBen bei jedem Rechnungs- 
schritt zu ermitteln. Einen brauchbaren Anhalt fiir die Konvergenz geben jedoch bereits die Langs- 
kraft oder das Biegemoment am Lochrand, die sich aus (42c) und (42d) leicht berechnen lassen, da 


n°) = N () und m\) ——— M¢) Hi Sane (43) 


Pn 
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Das Ergebnis einer solchen Rechnung ist in Tabelle 1 fiir fiinf verschiedene Werte &, angegeben. 
Dabei wird die gute Konvergenz des Verfahrens bestatigt: Bei kKleinen Léchern kann die Stérungs- 
rechnung somit bereits nach wenigen Schritten abgebrochen werden. 


Ne 


1 
S 
a 
x 


3,87 


— —— Soheibenlosung 


Lochrand €, 
Le 


; 


| 
| 


+5 


= YW 
> 


n) 
Abb.4. Verlauf der Funktionen - 7 (4 = 9,3). 


Tabelle 1. Langskrafte am Lochrand (u = 0,3). 
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e) Spannungs- und Formanderungszustand. Nach Abschlu8 der Stérungsrechnung 
lassen sich die SchnittgréBen aus (7) und (42) ermitteln. Mit 


co ive) 
== Se AG) = (») 
ae — ps i und Vo Ds, oe 
y=0 vy=0 


ergibt sich zum Beispiel die Langskraft 
a layie Ab L — “a ox as 
fee =2R BD [Amp pe + Bur Sn) + Bmp Zn) | (44) 
Pp 
wobei zur Abkirzung 


An,p = +2p(2p +1) = te sin ee + bn, p cos ae , 
Bap = 2 p(n 1) d,, p— (4n—4p +3) da p—1 + Gn, p=? > 
Bp = +2 (n—p) (2n—2p +1) dyp— (4n—4 p F 38) dy, pa + dnp 2 
gesetzt ist. 
Wahrend die Berechnung der zur Schalenmittelflache senkrechten Verschiebungen w unmittel- 
bar aus (42a) méglich ist, so ist die Ermittlung der beiden iibrigen Verschiebungskomponenten 
wesentlich umstandlicher. Zunachst scheint es naheliegend, w, aus (24a) durch Integration zu 


bestimmen. Dies ist jedoch unzweckmafig, da die dann auftretenden Integrale nicht in einfacher 
Form darstellbar sind. Es wird vielmehr aus (24b) und (24c) uw, eliminiert, und man erhalt 


Pst FRE IY ode eee 7 1 el cole 
Ora ey i &é Fern ae R (", rier gee net aie rns dé. (45) 
Dieses Integral kann geschlossen angegeben werden. So ist zum Beispiel fiir n = 0 


b= oR eB ae S he AG hee Z,.(€))} (46) 


Ist v,, ermittelt, so ]aBt sich wu, aus (24b) berechnen: 


(47) 


re é ns POD 1 1 1 
Uae “O & Eon a E Ca IR (", 9 Wr—l oe ay tonsa} 


SchlieBlich liefert (24a) eine Kontrolle fiir die Richtigkeit der Rechnung. 


6. Geschlossene Lésung fiir sehr kleine Lochdurchmesser. Im Bereich der Randstérung sei 


Gee Als (48a) 
nach (2) und (4) entspricht dies 
a< = ; (48b) 
Die Stérungsrechnung kann dann auf den ersten Rechnungsschritt beschrankt werden, also 
ye PO snp, (49) 


und die Hankelschen Funktionen lassen sich durch die ersten Glieder ihrer Potenzreihenentwicklung 


ersetzen. Nach einiger Zwischenrechnung erhalt man fiir die SchnittgréBen und Verschiebungen 
am Lochrand 


&2 
 &2 : 
Ny (¢=£,) & — (4 a. =| rsin2@, 


Lees g 
Me (Eg) & + FI ; z f (1 — p) In A Tsin29, 
oF me E2 - 
went) © +a (4 -+ Z| Tsin2, (50) 


Wet) = 1a 6p Ge sm 2 on 


Hierbei bezeichnet In y = 0,577216 die Eulersche Konstante. Ist §& = 0 und R = ov, so geht 
(50) in die Scheibenlésung iiber. Im iibrigen entsprechen die Gleichungen in ihrem Aufbau den 
von A. I. Lurje angegebenen Formeln fiir die auf axialen Zug beanspruchte Zylinderschale. 
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7. Vergleich der Spannungszustiinde bei Schale und Scheibe. Abb. 5 zeigt den Verlauf der 
Hauptspannung 
1 (m9— ng? o\ 


der Schalenmittelflache fiir 4 = 0,3 und den Lochhalbmesser & = 1,904, bezogen auf die Schub- 
spannung 


To = ie 
der ungeschwachten Schale. Abb. 6 
zeigt zum Vergleich die entspre- 
chende Darstellung fiir eine 
Scheibe. Dabei ist 


Q é a 
1) => ip mit No = “;. 


als dimensionslose Veranderliche 
eingefiihrt, wobei L eine beliebig 
wahlbare Lange, zum Beispiel die 
Langeneinheit, bezeichnet. 

Aus Abb.5 und 6 folgt, daB 
die Lochrandspannungen bei der 
Schale erheblich héher sind als bei 
der Scheibe. Wahrend die Rand- 
stérungen der Scheibe asympto- 
tisch gegen Null gehen, klingen sie 
bei der Schale in Form gedampfter 
Schwingungen ab. Die Abklingung 
ist dabei in Lochnahe wesentlich 
eréBer als bei der Scheibe. In wei- Abb. 5. Hauptspannung 
terer Entfernung vom Lochrand 
sind jedoch immer noch gering- 
fugige Spannungen vorhanden, die % | 
értlich sogar iiber den entspre- 
chenden Werten der Scheibe liegen 
kénnen. Diese schwache Damp- 
fung wurde von G. Sonntag! bei 
ahnlich belasteten Zylinderschalen 
ebenfalls festgestellt. 


oll 
qT 


der Schalenmittelfliche (4 = 0,3). 


KO 


Oa 


0 


8. Maximale Lochrandspannun- 
gen bei beliebigen Schalenabmes- 
sungen. Abb. 7 gibt den Verlauf 
der Normalkrafte n, am Lochrand 
fir die Halbmesser & < 4,0 an. 
Die Koordinate &;=0 entspricht 
der Scheibenlésung. Dabei ist 
beachtenswert, daB fiir gréBere &) 
die Stelle der maximalen Langs- ; Abb. 6. Hauptspannung 
kraft sich immer starker zur 
y-Achse hin verlagert. Die Spannungen werden also von der in ¢-Richtung starkeren Kriim- 
mung angezogen, was auch anschaulich unmittelbar einzusehen ist. 

Die Spannungsspitze in der Schalenleibung kann mit Hilfe von Abb. 8 unmittelbar angegeben 


werden, wobei der Einflu8 der Biegemomente beriicksichtigt wird. 


oT 


I 
der Scheibe. 


T) 


9. Zusammenfassung. Das Randstérungsproblem der tordierten Zylinderschale mit kreis- 
férmigem Loch wird mit Hilfe der Stérungsrechnung gelést. Die Lésungen konvergieren in dem 


1 G, Sonntag, Z. angew. Math, Mech. 36 (1956), S. 289. 
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fiir praktische Zwecke wesentlichen Parameterbereich sehr gut. Es wird das ae aeaeher 
gefundene Ergebnis bestatigt, wonach die Randspannungen bei der Schale stets gré : ae als 
bei der Scheibe. Diese Spannungserhéhung rithrt nicht nur von den zusitzlich auftreten en Biege- 
momenten her, sondern beruht auch auf einem Anwachsen der Langskrafte. Die Stérungen papi 
in Form gedaimpfter Schwingungen ab, die sich mathematisch durch negative Potenzreihen un 


Gz 


O0ey 


! 
0 005“ O07 Qe 005 07 02 04 


Onax a 26 1% 


Abb. 8. Spannungserhéhung am Lochrand (i = 0,3). 


Reihen Hankelscher Zylinderfunktionen ausdriicken lassen. Die SchnittgréBen und Verformungen 
des Lochrandes werden fiir sehr kleine Lochdurchmesser in geschlossener Form angegeben. Fiir 
beliebige Schalenabmessungen kann die Spannungserhéhung am Lochrand in einem Diagramm 
abgelesen werden. 


(Eingegangen am 8. Marz 1958.) 
Anschrift des Verfassers: Dr.-Ing. Dieter Withum, Duisburg-Huckingen, Am Heidberg 27. 
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